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Abstract 

A priori, there is nothing very special about shear-free or asymptotically shear-free null 
geodesic congruences. Surprisingly, however, they turn out to possess a large number of fasci- 
nating geometric properties and to be closely related, in the context of general relativity, to a 
variety of physically significant effects. It is the purpose of this paper to try to fully develop 
these issues. 

This work starts with a detailed exposition of the theory of shear-free and asymptotically 
shear-free null geodesic congruences, i.e., congruences with shear that vanishes at future con- 
formal null infinity. A major portion of the exposition lies in the analysis of the space of regular 
shear-free and asymptotically shear-free null geodesic congruences. This analysis leads to the 
space of complex analytic curves in an auxiliary four-complex dimensional space, "H-space. 
They in turn play a dominant role in the applications. 

The applications center around the problem of extracting interior physical properties of an 
asymptotically-flat spacetime directly from the asymptotic gravitational (and Maxwell) field 
itself, in analogy with the determination of total charge by an integral over the Maxwell field 
at infinity or the identification of the interior mass (and its loss) by (Bondi's) integrals of the 
Weyl tensor, also at infinity. 

More specifically, we will see that the asymptotically shear-free congruences lead us to 
an asymptotic definition of the center-of-mass and its equations of motion. This includes a 
kinematic meaning, in terms of the center-of-mass motion, for the Bondi three-momentum. In 
addition, we obtain insights into intrinsic spin and, in general, angular momentum, including 
an angular-momentum-conservation law with well-defined flux terms. When a Maxwell field 
is present, the asymptotically shear-free congruences allow us to determine/define at infinity a 
center-of-charge world line and intrinsic magnetic dipole moment. 
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1 Introduction 

Though from the very earliest days of Lorentzian geometries, families of null geodesies (null geodesic 
congruences (NGCs)) were obviously known to exist, it nevertheless took many years for their 
significance to be realized. It was from the seminal work of Bondi [8J, with the introduction 
of null surfaces and their associated null geodesies used for the study of gravitational radiation, 
that the realization of the importance of NGCs became recognized. To analyze the differential 
structure of such congruences, Sachs [56] introduced the fundamental 'tools', known as the optical 
parameters, namely, the divergence, the twist (or curl) and the shear of the congruence. From 
the optical parameters one then could classify congruences by the vanishing (or the asymptotic 
vanishing) of one or more of these parameters. All the different classes exist in flat space but, in 
general, only special classes exist in arbitrary spacetimes. For example, in flat space, divergence-free 
congruences always exist, but for nonflat vacuum spacetimes they exist only in the case of certain 
high symmetries. On the other hand, twist-free congruences (null surface-forming congruences) 
exist in all Lorentzian spacetimes. General vacuum spacetimes do not allow shear-free congruences, 
though all asymptotically-flat spacetimes do allow asymptotically shear-free congruences, a natural 
generalization of shear-free congruences, to exist. 

Our primary topic of study will be the cases of shear-free and asymptotically shear-free NGCs. In 
flat space the general shear-free congruences have been extensively studied. However, only recently 
has the special family of regular congruences been investigated. In general, as mentioned above, 
vacuum (or Einstein-Maxwell) metrics do not possess shear- free congruences; the exceptions being 
the algebraically-special metrics, all of which contain one or two such congruences. On the other 
hand, all asymptotically-flat spacetimes possess large numbers of regular asymptotically shear-free 
congruences. 

A priori there does not appear to be anything very special about shear-free or asymptotically 
shear-free NGCs. However, over the years, simply by observing a variety of topics, such as the clas- 
sification of Maxwell and gravitational fields (algebraically-special metrics), twistor theory, H-space 
theory and asymptotically-flat spacetimes, there have been more and more reasons to consider them 
to be of considerable importance. One of the earliest examples of this is Robinson's [53] demon- 
stration that a necessary condition for a curved spacetime to admit a null solution of Maxwell's 
equation is that there be, in that space, a congruence of null, shear-free geodesies. Recent results 
have shown that the regular congruences - both the shear-free and the asymptotically shear-free 
congruences - have certain very attractive and surprising properties; each congruence is determined 
by a complex analytic curve in the auxiliary complex space that is referred to as "H-space. For 
asymptotically-flat spacetimes, some of these curves contain a great deal of physical information 
about the spacetime itself [29l [27l [28] . 

It is the main purpose of this work to give a relatively complete discussion of these issues. 
However, to do so requires a digression. 

A major research topic in general relativity (GR) for many years has been the study of asymptotically- 
flat spacetimes. Originally, the term 'asymptotically flat' was associated with gravitational fields, 
arising from finite bounded sources, where infinity was approached along spacelike directions (e.g., [5] 
[58]). Then the very beautiful work of Bondi [8 showed that a richer and more meaningful idea to 
be associated with 'asymptotically flat' was to study gravitational fields in which infinity was ap- 
proached along null directions. This led to an understanding of gravitational radiation via the Bondi 
energy-momentum loss theorem, one of the profound results in GR. The Bondi energy-momentum 
loss theorem, in turn, was the catalyst for the entire contemporary subject of gravitational radia- 
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tion and gravitational wave detectors. The fuzzy idea of where and what is infinity was clarified 
and made more specific by the work of Penrose [46J [47] with the introduction of the conformal 
compactification (via the rescaling of the metric) of spacetime, whereby infinity was added as a 
boundary and brought into a finite spacetime region. Penrose's infinity or spacetime boundary, 
referred to as Scri or 3, has many sub-regions: future null infinity, 3 + ; past null infinity, 3~; future 
and past timelike infinity, I + and I ; and spacelike infinity, 7°. In the present work, 3 + and its 
neighborhood will be our arena for study. 

A basic question for us is what information about the interior of the spacetime can be obtained 
from a study of the asymptotic gravitational field; that is, what can be learned from the remnant 
of the full field that now 'lives' or is determined on 3 + ? This quest is analogous to obtaining the 
total interior electric charge or the electromagnetic multipolc moments directly from the asymptotic 
Maxwell field, i.e., the Maxwell field at 3 + , or the Bondi energy- momentum four- vector from the 
gravitational field (Weyl tensor) at 3 + . However, the ideas described and developed here are not 
in the mainstream of GR; they lie outside the usual interest and knowledge of most GR workers. 
Nevertheless, they are strictly within GR; no new physics is introduced; only the vacuum Einstein 
or Einstein-Maxwell equations are used. The ideas come simply from observing (discovering) cer- 
tain unusual and previously overlooked features of solutions to the Einstein equations and their 
asymptotic behavior. 

These observations, as mentioned earlier, centered on the awakening realization of the remarkable 
properties and importance of the special families of null geodesies: the regular shear-free and 
asymptotically shear-free NGCs. 

The most crucial and striking of these overlooked features (mentioned now but fully developed 
later) are the following: in flat space every regular shear-free NGC is determined by the arbitrary 
choice of a complex analytic world line in complex Minkowski space, Mc< Furthermore and more 
surprising, for every asymptotically-flat spacetime, every regular asymptotically shear-free NGC 
is determined by the given Bondi shear (given for the spacetime itself) and by the choice of an 
arbitrary complex analytic world line in an auxiliary complex four-dimensional space, 'H-space, 
endowed with a complex Ricci-flat metric structure. In other words, the space of regular shear-free 
and asymptotically shear-free NGCs are both determined by arbitrary analytic curves in Mc and 
H-space respectively [29 | [27 ] [26] . 

Eventually, a unique complex world line in this space is singled out, with both the real and 
imaginary parts being given physical meaning. The detailed explanation for the determination of 
this world line is technical and reserved for a later discussion. However, a rough intuitive idea can 
be given in the following manner. 

The idea is a generalization of the trivial procedure in electrostatics of first defining the electric 
dipole moment, relative to an origin, and then shifting the origin so that the dipole moment vanishes 
and thus obtaining the center of charge. Instead, we define, on 3 + , with specific Bondi coordinates 
and tetrad, the complex mass dipole moment (the real mass dipole plus 'i' times angular momentum) 
from certain components of the asymptotic Weyl tensor. (The choice of the specific Bondi system 
is the analogue of the choice of origin in the electrostatic case.) Then, knowing how the asymptotic 
Weyl tensor transforms under a change of tetrad and coordinates, one sees how the complex mass 
dipole moment changes when the tetrad is rotated to one defined from the asymptotically shear-free 
congruence. By setting the transformed complex mass dipole moment to zero, the unique complex 
world line, identified as the complex center of mass, is obtained. (In Einstein-Maxwell theory a 
similar thing is done with the asymptotic Maxwell field leading to the vanishing of the complex 
Maxwell dipole moment [electric plus '«' times magnetic dipole moment], with a resulting complex 
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center of charge.) 

This procedure, certainly unusual and out of the mainstream and perhaps appearing to be 
ambiguous, does logically hold together. The real justification for these identifications comes not 
from its logical structure, but rather from the observed equivalence of the derived results from these 
identifications with well-known classical mechanical and electrodynamical relations. These derived 
results involve both kinematical and dynamical relations. Though they will be discussed at length 
later, we mention that they range from a kinematic expression for the Bondi momentum of the 
form, P = Mv + . . .; a derivation of Newton's second law, F = Ma; a conservation law of angular 
momentum with a well-known angular momentum flux; to the prediction of the Dirac value of the 
gyromagnetic ratio. We note that, for the charged spinning particle metric [37], the imaginary part 
of the world line is indeed the spin angular momentum, a special case of our results. 

A major early clue that shear-free NGCs were important in GR was the discovery of the (vacuum 
or Einstein-Maxwell) algebraically special metrics. These metrics are defined by the algebraic 
degeneracy in their principle null vectors, which form (by the Goldberg-Sachs theorem [18]) a null 
congruence which is both geodesic and shear-free. For the asymptotically-flat algebraically-special 
metrics, this shear- free congruence (a very special congruence from the set of asymptotically shear- 
free congruences) determines a unique world line in the associated auxiliary complex "H-space. This 
shear-free congruence (with its associated complex world line) is a special case of the above argument 
of transforming to the complex center of mass. Our general asymptotically-flat situation is, thus, a 
generalization of the algebraically-special case. Much of the analysis leading to the transformation 
of the complex dipoles in the case of the general asymptotically-flat spaces arose from generalizing 
the case of the algebraically-special metrics. 

To get a rough feeling (first in flat space) of how the curves in Mc are connected with the 
shear-free congruences, we first point out that the shear-free congruences are split into two classes: 
the twisting congruences and the twist-free ones. The regular twist-free ones are simply the null 
geodesies (the generators) of the light cones with apex on an arbitrary timelike Minkowski space 
world line. Observing backwards along these geodesies from afar, one 'sees' the world line. The 
regular twisting congruences are generated in the following manner: consider the complexification 
of Minkowski space, Mc- Choose an arbitrary complex (analytic) world line in Mc and construct 
its family of complex light cones. The projection into the real Minkowski space, M, of the complex 
geodesies (the generators of these complex cones), yields the real shear-free twisting NGCs [4]. The 
twist contains or 'remembers' the apex on the complex world line. Looking backwards via these 
geodesies, one appears 'to see' the complex world line. In the case of asymptotically shear-free 
congruences in curved spacetimes, one can not trace the geodesies back to a complex world line. 
However, one can have the illusion (i.e., a virtual image) that the congruence is coming from a 
complex world line. It is from this property that we can refer to the asymptotically shear-free 
congruences as lying on generalized light cones. 

The analysis of the geometry of the asymptotically shear-free NGCs is greatly facilitated by the 
introduction of Good-Cut Functions (GCFs). Each GCF is a complex slicing of 3 + from which the 
associated asymptotically shear-free NGC and world line can be easily obtained. For the special 
world line and congruence that leads to the complex center of mass, there is a unique GCF that is 
referred to as the Universal-Cut Function, (UCF). 

Information about a variety of objects is contained in and can be easily calculated from the UCF, 
e.g., the unique complex world line; the direction of each geodesic of the congruence arriving at 
3+; and the Bondi asymptotic shear of the spacetime. The ideas behind the GCFs and UCF arose 
from some very pretty mathematics: from the 'good-cut equation' and its complex four-dimensional 
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solution space, "H-space [331 123] ■ In flat space almost every asymptotically vanishing Maxwell field 
determines its own Universal Cut Function, where the associated world line determines both the 
center of charge and the magnetic dipole moment. In general, for Einstein-Maxwell fields, there 
will be two different UCFs, (and hence two different world lines), one for the Maxwell field and one 
for the gravitational field. The very physically interesting special case where the two world lines 
coincide will be discussed. 

In this work, we seek to provide a comprehensive overview of the theory of asymptotically shear- 
free NGCs, as well as their physical applications to both flat and asymptotically-flat spacetimes. 
The resulting theoretical framework unites ideas from many areas of relativistic physics and has 
a crossover with several areas of mathematics, which had previously appeared short of physical 
applications. 

The main mathematical tool used in our description of 3 + is the Newman-Penrose (NP) formal- 
ism |39j . Spherical functions are expanded in spin-s tensor harmonics |43) ; in our approximations 
only the I = 0, 1, 2 harmonics are retained. Basically, the detailed calculations should be considered 
as expansions around the Reissner-Nordstrom metric, which is treated as zeroth order; all other 
terms being small, i.e., at least first order. We retain terms only to second order. 

In Section [3J we give a brief review of Penrose's conformal null infinity 3 along with an expo- 
sition of the NP formalism and its application to asymptotically-flat spacetimes. There is then a 
description of J + , the stage on which most of our calculations take place. The Bondi mass aspect 
(a function on 3 + ) is defined from the asymptotic Weyl tensor. From it we obtain the physical 
identifications of the Bondi mass and linear momentum. Also discussed is the asymptotic sym- 
metry group of 3 + , the Bondi-Metzner-Sachs (BMS) group [51 [551 [301 [3H|- The Bondi mass and 
linear momentum become basic for the physical identification of the complex center-of-mass world 
line. For its pedagogical value and prominence in what follows, we review Maxwell theory in the 
spin-coefficient (SC) formalism. 

Section [3] contains the detailed analysis of shear-free NGCs in Minkowski spacetime. This 
includes the identification of the flat space GCFs from which all regular shear-free congruences 
can be found. We also show the intimate connection between the flat space GCFs, the (homo- 
geneous) good-cut equation, and Mc- As applications, we investigate the UCF associated with 
asymptotically-vanishing Maxwell fields and in particular the shear-free congruences associated 
with the Lienard-Wiechert (and complex Lienard-Wiechert) fields. This allows us to identify a real 
(and complex) center-of-charge world line, as mentioned earlier. 

In Section^ we give an overview of the machinery necessary to deal with twisting asymptotically 
shear-free NGCs in asymptotically-flat spacetimes. This involves a discussion of the theory of H- 
space, the construction of the good-cut equation from the asymptotic Bondi shear and its complex 
four-parameter family of solutions. We point out how the simple Minkowski space of the preceding 
Section [3] can be seen as a special case of the more general theory outlined here. These results have 
ties to Penrose's twistor theory and the theory of Cauchy-Riemann (CR) structures; an explanation 
of these crossovers is given in Appendicies IAIB1 

In Section [5j as examples of the ideas developed here, linear perturbations off the Schwarz- 
schild metric, the algebraically-special type II metrics and asymptotically-stationary spacetime are 
discussed. 

In Section [51 the ideas laid out in the previous Sections [3J [3] and [5J are applied to the general 
class of asymptotically-flat spacetimes: vacuum and Einstein-Maxwell. Here, reviewing the material 
of the previous section, we apply the solutions to the good-cut equation to determine all regular 
asymptotically shear-free NGCs by first choosing arbitrary world lines in the solution space and then 
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singling out a unique one; two world lines in the Einstein-Maxwell case, one for the gravitational 
field, the other for the Maxwell field. This identification of the unique lines comes from a study of 
the transformation properties, at 3 + , of the asymptotically-defined mass and spin dipoles and the 
electric and magnetic dipoles. The work of Bondi, with the identification of energy-momentum and 
its evolution, allows us to make a series of surprising further physical identifications and predictions. 
In addition, with a slightly different approximation scheme, we discuss our ideas applied to the 
asymptotic gravitational field with an electromagnetic dipole field as the source. 

Section [7] contains an analysis of the gauge (or BMS) invariance of our results. 

Sectional the Discussion/ Conclusion section, begins with a brief history of the origin of the ideas 
developed here, followed by comments on alternative approaches, possible physical predictions from 
our results, a summary and open questions. 

Finally, we conclude with three appendices, which contain several mathematical crossovers that 
were frequently used or referred to in the text: CR structures and twistors, along with a brief expo- 
sition of the tensorial spherical harmonics |43j and their Clebsch-Gordon product decompositions. 

1.1 Notation and definitions 

The following contains the notational conventions that will be in use throughout the course of this 
review. 

• Throughout this work we use the symbols T, 'to', 'n' . . .with several different 'decorations' 
but always meaning a null tetrad or a null tetrad field. 

a) Though in places, e.g., in Section [2~4l the symbols, l a , m a , n a . . . , i.e., with an a, 6,c. . . 
can be thought of as the abstract representation of a null tetrad (i.e., Penrose's abstract index 
notation [50] ) . in general, our intention is to describe vectors in a coordinate representation. 

b) The symbols, l a , l# a ,l* a most often represent the coordinate versions of different null 
geodesic tangent fields, e.g., one- leg of a Bondi tetrad field or some rotated version. 

c) The symbol, l a , (with hat) has a very different meaning from the others. It is used to 
represent the Minkowski components of a normalized null vector giving the null directions on 
an arbitrary light cone: 

v2 (i+cc,c+c,ic-ic,-i+cc) = f^o >k°V a) 



2(1 + CO V 2 2 

As the complex stereographic coordinates (C, sweep out the sphere, the l a sweeps out the 
entire future null cone. The other members of the associated null tetrad are 



V2 



2(1 + CO 
2(1 + CO 



( ,l-C 2 ,- l (l + C 2 ),2C), (2) 

(i + CC,-(C + OX-*C,i-CC)- 



• Several different timelike variables (mb, u rc t, T, s) and their derivatives are used. 
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The Bondi time, ub, is closely related to the retarded time, u re t = \/2ub- The use of the 
retarded time, u re t> is important in order to obtain the correct numerical factors in the ex- 
pressions for the final physical results. Their derivatives are represented by 

d UB K = K, (3) 

d Utet K = K' = ^K. 

The u rc t,T,s, derivatives are denoted by the same prime (') since it is always applied to 
functions with the same functional argument. Though we are interested in real physical 
spacetime, often the time variables (u TCt , ub, t) take complex values close to the real (s is 
always real). Rather than putting on 'decorations' to indicate when they are real or complex 
(which burdens the expressions with an overabundance of different symbols), we leave reality 
decisions to be understood from context. In a few places where the reality of the particular 
variable is manifestly first introduced (and is basic) we decorate the symbol by a superscript 
(i?), i.e., or t^et ■ After their introduction we revert to the undecorated symbol. 

Note: At this point we are taking the velocity of light as c = 1 and omitting it; later, when 
we want the correct units to appear explicitly, we restore the c. This entails, via r —¥ or, 
s — > cs, changing the prime derivatives to include the c, i.e., 

K' -> c^K'. (4) 

• Often the angular (or sphere) derivatives, S and 3, are used. The notation 3( a )-fC means, 
apply the S operator to the function K while holding the variable (a) constant. 

• The complex conjugate is represented by the overbar, e.g., £. When a complex variable, (, is 
close to the complex conjugate of £, but independent, we use £ s» (. 

Frequently, in this work, we use terms that are not in standard use. It seems useful for clarity 
to have some of these terms defined early. 

• As mentioned earlier, we use the term "generalized light cones" to mean (real) NGCs that 
appear to have their apexes on a world line in the complexification of the spacetime. A detailed 
discussion of this will be given in Sections [3] and |U 

• The term "complex center of mass" (or "complex center of charge") is frequently used. Up 
to the choice of constants (to give correct units) this is basically the "real center of mass plus 

angular momentum" (or "real center of charge plus H' magnetic moment"). There will 
be two different types of these "complex centers . . . " ; one will be geometrically defined or 
intrinsic, i.e., independent of the choice of coordinate system, the other will be relative, i.e., it 
will depend on the choice of (Bondi) coordinates. The relations between them are nonlinear 
and nonlocal. 

• A very important technical tool used throughout this work is a class of complex analytic 
functions, ub — G(t, £, £), referred to as GCFs that are closely associated with shear-free 
NGCs. The details are given later. For any given asymptotically-vanishing Maxwell field with 
nonvanishing total charge, the Maxwell field itself allows one, on physical grounds, to choose a 
unique member of the class referred to as the (Maxwell) UCF. For vacuum asymptotically-flat 
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spacetimes, the Weyl tensor allows the choice of a unique member of the class referred to as 
the (gravitational) UCF. For Einstein-Maxwell there will be two such functions, though in 
important cases they will coincide and be referred to as UCFs. When there is no ambiguity, 
in either case, they will simple be UCFs. 

• A notational irritant arises from the following situation. Very often we expand functions on 
the sphere in spin-s harmonics, as, e.g., 

X = X°Y a + x'YuiCO + X ij Y 2ij ((,0 + X ijk Y 3ij k((,() + ■■-, 

where the indices, i,j,k... represent three-dimensional Euclidean indices. To avoid extra 
notation and symbols we write scalar products and cross-products without the use of an 
explicit Euclidean metric, leading to awkward expressions like 

if = rfX^rfXi, 

H k = (if x X) k sfj'A^eyfc. 

This, though easy to understand and keep track of, does run into the unpleasant fact that 
often the relativist four-vector, 

x (x 7 X )? 

appears as the I = 0, 1 harmonics in the harmonic expansions. Thus, care must be used when 
lowering or raising the relativistic index, i.e., r\ a bX a = Xb = (x°> ~ X*)- 

We conclude this introduction with a comprehensive glossary of symbols and terms used throughout 
the text. 
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Tabic 1: Glossary 



Symbol/ Acronym 



M, M c 

MB, "rot 
Sub / = / 

r 

(CO 

V?J..,y(C-C) 



9(a)/ 

{/ a ,n a ,m a ,m a } 

NGC 

NP/SC 

{lp ,1pl,1p2,1pS,1p4} 
{00,01,02} 

P 

s 

k 

u = g(t , c, C) 

r = s + iX = T(u, C, C) 

d T f = f 

3f T) G(r,C,C) = ^°(«.C,0 
GCF 

L( u ,C,C)=3(r)G 

g (MB) T + LT = 
three-dimensional CR Structure 
"H-space 



Definition 



Future conformal null infinity, Complex future conformal null 
infinity 

Future, Past conformal timelike infinity, Conformal spacelike 
infinity 

Minkowski space, Complex Minkowski space 

Bondi time coordinate, Retarded Bondi time (\/2ub = "rot) 

Derivation with respect to ub 

Derivation with respect to w ro t 

Affine parameter along null geodesies 

(e^cot(0/2), e _i * cot(0/2)); stereographic coordinates on S 2 
Tensorial spin-s spherical harmonics 

P 1s ~§qP s -> P 1+s -§qP~ s \ spin-weighted operator on the two- 
sphere 

Metric function on S 2 ; often P = P = 1 + (( 

Application of 5-operator to / while the variable a is held 

constant 

Null tetrad system; l a n a = —m a fh a = 1 
Null Geodesic Congruence 
Newman-Penrose/Spin-Coefficient Formalism 
Metric coefficients in the Newman-Penrose formalism 
Weyl tensor components in the Newman-Penrose formalism 
Maxwell tensor components in the Newman-Penrose formal- 
ism 

Complex divergence of a null geodesic congruence 

Twist of a null geodesic congruence 

Complex shear, Asymptotic complex shear of a NGC 

^F-; Gravitational constant 

Cut function on 3 + 

Complex auxiliary (CR) potential function 
Derivation with respect to r 

Good-Cut Equation, describing asymptotically shear-free 
NGCs 

Good-Cut Function 

Stereographic angle field for an asymptotically shear-free NGC 
at 3+ 

CR equation, describing the embedding of 3 + into C 2 

A class of one-forms describing a real three-manifold of C 2 

Complex four-dimensional solution space to the Good-Cut 

Equation 
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Symbol/ Acronym 



Definition 



^C(grav) - 



iD 



M 



(mass) 1 



UCF 

u = X(t,C,0 



* = ^ - 

M B = - 
F 



S 2 (J + (7(7 
$0 



2V2G 

- e 9 



Complex electromagnetic dipole 

Complex center-of-charge world line, lives in "H-space 
Complex gravitational dipole 

Complex center-of-mass world line, lives in H-space 
Universal Cut Function 

UCF; corresponding to the complex center-of-charge world line 
Bondi Mass Aspect 
Bondi mass 

Bondi linear three-momentum 

Vacuum linear theory identification of angular momentum 
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2 Foundations 

In this section, we review several of the key ideas and tools that are indispensable in our later 
discussions. We keep our explanations as concise as possible, and refrain from extensive proofs of 
any propositions. The reader will be directed to the appropriate references for the details. In large 
part, much of what is covered in this section should be familiar to many workers in GR. 

2.1 Asymptotic flatness and J + 

Ever since the work of Bondi [5] illustrated the importance of null hypersurfaces in the study of 
outgoing gravitational radiation, the study of asymptotically-flat spacetimes has been one of the 
more important research topics in GR. Qualitatively speaking, a spacetime can be thought of as 
(future) asymptotically flat if the curvature tensor vanishes at an appropriate rate as infinity is 
approached along the future-directed null geodesies of the null hypersurfaces. The type of physical 
situation we have in mind is an arbitrary compact gravitating source (perhaps with an electric 
charge and current distribution), with the associated gravitational (and electromagnetic) field. The 
task is to gain information about the interior of the spacetime from the study of far-field features, 
multipole moments, gravitational and electromagnetic radiation, etc. |44j . The arena for this study 
is on what is referred to as future null infinity, 3 + , the future boundary of the spacetime. The 
intuitive picture of this boundary is the set of all endpoints of future-directed null geodesies. 

A precise definition of null asymptotic flatness and the boundary was given by Penrose [JHl HZ] j 
whose basic idea was to rescale the spacetime metric by a conformal factor, which approaches zero 
asymptotically: the zero value defining future null infinity. This process leads to the boundary being 
a null hypersurface for the conformally-rescaled metric. When this boundary can be attached to the 
interior of the rescaled manifold in a regular way, then the spacetime is said to be asymptotically 
flat. 

As the details of this formal structure are not used here, we will rely largely on the intuitive 
picture. A thorough review of this subject can be found in |12) . However, there are a number of 
important properties of 3 + arising from Penrose's construction that we rely on [44j |46j [47] : 

(A) : For both the asymptotically- flat vacuum Einstein equations and the Einstein-Maxwell 
equations, J + is a null hypersurface of the conformally rescaled metric. 

(B) : 3+ is topological^ S 2 x M. 

(C) : The Weyl tensor C£ cd vanishes at 3 + , with the peeling theorem describing the speed of its 
falloff (see below). 

Property (B) allows an easy visualization of the boundary, 3 + , as the past light cone of the 
point J + , future timelike infinity. As mentioned earlier, J + will be the stage for our study of 
asymptotically shear-free NGCs. 

2.2 Bondi coordinates and null tetrad 

Proceeding with our examination of the properties of 3 + , we introduce, in the neighborhood of 
3 + , what is known as a Bondi coordinate system: {ub,t, C>C)- I n this system, ub, the Bondi 
time, labels the null surfaces, r is the affine parameter along the null geodesies of the constant ub 
surfaces and £ = e 1 ^ cot(0/2), the complex stereographic angle labeling the null geodesies of 3 + . 
To reach J+ we simply let r — > oo, so that 3 + has coordinates (ub,C>0- The time coordinate 
mb, the topologically M portion of 3 + , labels "cuts" of 3 + . The stereographic angle ( accounts for 



2 FOUNDATIONS 



14 



the topological generators of the S 2 portion of 3 + , i.e., the null generators of 3 + . The choice of a 
Bondi coordinate system is not unique, there being a variety of Bondi coordinate systems to choose 
from. The coordinate transformations between any two, known as Bondi-Metzner-Sachs (BMS) 
transformations or as the BMS group, are discussed later in this section. 

Associated with the Bondi coordinates is a (Bondi) null tetrad system, (l a , n a , m a ,m a ) . The 
first tetrad vector l a is the tangent to the geodesies of the constant uq null surfaces given by [4"4"] 



V 

l a VJ b 

d 



l a 



dx a 



dx a 
dr 

0. 
d_ 

dr 



9 ab V b u B , 



The second null vector n a is normalized so that: 
In Bondi coordinates, we have [44] 



1. 



d 



8 



dx a dun 



U 



di 



-X J 



dx A 



(5) 
(6) 
(7) 

(8) 
(9) 



for functions U and X A to be determined, and A — 2,3 with x 2 = £ and a; 3 = C- At 3 + , n a is 
tangent to the null generators of J + . 

The tetrad is completed with the choice of a complex null vector m a , (m a m a — 0) which is 
itself orthogonal to both l a and n ai initially tangent to the constant ub cuts at 3 + and parallel 
propagated inward on the null geodesies. It is normalized by 



771 771 o 



-1. 



(10) 



Once more, in coordinates, we have 1441 



m 



dx a 



di 



dx A ' 



for some uj and £ A to be determined. All other scalar products in the tetrad are to vanish. 
With the tetrad thus defined, the contravariant metric of the spacetime is given by 



-itb 



g~" = l a n° + l b n a - m~m" 
In terms of the metric coefficients U, uj, X a , and £ A , the metric can be written as: 



m b m a 



9 = 



1 





l A 



AB 





1 9" 9 
g 1A 9 AB 

2(U-uQ), 



(11) 



(12) 



(13) 
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Thus, the spacetime metric is determined from the metric coefficients. 

There remains the issue of both coordinate and tetrad freedom, i.e., local Lorentz transforma- 
tions. Most of the time we work in one arbitrary but fixed Bondi coordinate system, though for 
special situations more general coordinate systems are used. The more general transformations 
are given, essentially, by choosing an arbitrary slicing of 3 + , written as ub = G(s, C, C) with s 
labeling the slices. To keep conventional coordinate conditions unchanged requires a rescaling of 
r : r — s- r' = (d s G)~ 1 r. It is also useful to be able to shift the origin of r by r' = r — ro(ws, C> C) 
with arbitrary to(ub, C> C)- 

The tetrad freedom of null rotations around n a , performed in the neighborhood of 3 + , will later 
play a major role. For an arbitrary function L(ub, C, C) on 3 + , the null rotation about the vector 
n a [33] is given by 

l a -> r a = r--m Q --m a + 0(r~ 2 ), (14) 

r 





L 




-> r a = i a - 


r 


m a 


— > m* a = m a 




L 




—'i 
r 








-> n* a =n a . 







m a -> m* a = m a - -n a + 0(r" 2 ), (15) 
r 

n a -> n* a =n a . (16) 

Eventually, by the appropriate choice of the function L(ub,CiC)j the- new null vector, l* a , can 
be made into the tangent vector of an asymptotically shear-free NGC. 

The second type of tetrad transformation is the rotation in the tangent (m a , m a ) plane, which 
keeps l a and n a fixed: 

m a -> e iX m a , A e M. (17) 

This latter transformation provides motivation for the concept of spin weight. A quantity r]r s -\ (£, Q 
is said to have spin- weight s if, under the transformation, Equation (I17|) . it transforms as 

V^VU(C,<) =e*% (£C)- (18) 

An example would be to take a vector on 3 + , say rj a , and form the spin- weight-one quantity, 

V(i) = ^ma- 
Comment: For later use we note that L(ub, Cj C) has spin weight, s = 1. 

For each s, spin-s functions can be expanded in a complete basis set, the spin-s harmonics, 

sYi m (C,C) or spin-s tensor harmonics, ^/^(C »*5m(C>C)- 
A third tetrad transformation, the boosts, are given by 

l* a = Kl c \ n* a = K- Y n a . (19) 

These transformations induce the idea of conformal weight, an idea similar to spin weight. Under 
a boost transformation, a quantity, r\i w ) , will have conformal weight w if 

V^)^vt ) =K v, V(v,y (20) 

Sphere derivatives of spin- weighted functions rj^ (£, C) are given by the action of the operators 
8 and its conjugate operator 5, defined by [17] 

^=P 1 -^^, (2D 
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%)=P 1+S ^^, (22) 



where the function P is the conformal factor defining the sphere metric, 

ds 2 



2 _ 



P 2 

most often taken as 



p = p = i + CC- 



2.3 The optical equations 



Since this work concerns NGCs and, in particular, shear-free and asymptotically shear-free NGCs, 
it is necessary to first define them and then study their properties. 

Given a Lorentzian manifold with local coordinates, x a , and an NGC, i.e., a foliation by a three 
parameter family of null geodesies, 

x a = X a (r,y w ), (23) 

with r the affine parameterization and the (three) y w labeling the geodesies, the tangent vector field 
l a = DX a = d r X a satisfies the geodesic equation 

l a V a l b = 0. 

The two complex optical scalars (spin coefficients), p and a, arc defined by 

P = \(-VJ a +i curlH, (24) 
curU a = y/(V [a l b ]V a l b ) 

and 

o- = V [a l b) m a m b 

with m a an arbitrary complex (spacelike) vector satisfying m a m a = m a l a = m a m a — 1 = 0. 
Equivalently a can be defined by its norm, 



1 / _ . _.. ., J 

2 



^=dV ( A)V^--(V/) 2 



with an arbitrary phase. 

The p and a satisfy the optical equations of Sachs [56], namely, 

dp _ 2 



p 2 + aa + $ 00 , (25) 
-^- = 2p<j + ^ (26) 

^00 = Rab l a l b , 

= -C a bc d l a m b l c m d , 
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where $00 and 4"o are, respectively, a Ricci and a Weyl tensor tetrad component (see below). In flat 
space, i.e., with $00 = ^0 = 0, excluding the degenerate case of pp — 00 = 0, plane and cylindrical 
fronts, the general solution is 

zE — r 

P= T* + W-oW (27) 

5r (28) 



a 



r 2 + E 2 - crV 

The complex cr° (referred to as the asymptotic shear) and the real S (called the twist) are determined 
from the original congruence, Equation (|2"5|) . Both are functions just of the parameters, y w . Their 
behavior for large r is given by 

1 iE E 2 a°a° 



0(r- 4 ), (29) 



o = ^ + 0(r- 4 ). (30) 

From this, er° gets its name as the asymptotic shear. In Section [31 we return to the issue of the 
explicit construction of NGCs in Minkowski space and in particular to the construction and detailed 
properties of regular shear-free congruences. 

Note the important point that, in M, the vanishing of the asymptotic shear forces the shear 
to vanish. The same is not true for asymptotically-flat spacetimes. Specifically, for future null 
asymptotically-flat spaces described in a Bondi tetrad and coordinate system, we have, from other 
considerations, that 

<foo = 0(r~ e ), 
^0 = 0(r~ 5 ), 
E = 0, 

which leads to the asymptotic behavior of p and er, 

p = P=-- + — + 0(r ), 
o = °f 2+ 0{r-% 

with the two order symbols explicitly depending on the leading terms in $00 and ^0. The vanishing 
of do does not, in this nonflat case, imply that a vanishes. This case, referred to as asymptotically 
shear-free, plays the major role later. It will be returned to in greater detail in Section 2] 



2.4 The Newman— Penrose formalism 

Though the NP formalism is the basic working tool for our analysis, this is not the appropriate 
venue for its detailed exposition. Instead we will simply give an outline of the basic ideas followed 
by the results found, from the application of the NP equations, to the problem of asymptotically-flat 
spacetimes. 

The NP version [39l [44j [41] of the vacuum Einstein (or the Einstein-Maxwell) equations uses 
the tetrad components 

\? = (l a ,n a ,m a ,rh a ), (31) 
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(i = 1, 2, 3, 4) rather than the metric, as the basic variable. (An alternate version, not discussed here, 
is to use a pair of two-component spinors.) The metric, Equation (|12[) . can be written compactly 
as 



with 





(32) 


/ 1 \ 




10 


(33) 


-1 


\ -1 J 





The complex spin coefficients, which play the role of the connection, are determined from the 
Ricci rotation coefficients [39l [44] : 

l)k = A j A ffla A l = A jK;k> (34) 

via the linear combinations 

a = | (7124 - 7344), A 
^=9(7123-7343), fJ- 
7 = 9(7122 - 7342), v 

£=2(7l21 - 734l), 7T 

The third basic variable in the NP formalism is the Weyl tensor or, equivalently, the following 
five complex tetrad components of the Weyl tensor: 

4>q = -C abcd l a m b l c m d , fa = -C abcd l a n b l c m d , (36) 

i'2 = -\ (Cabcd l a n b l c n d - C abcd l a n b m c m d ) , (37) 

V-3 = C abcd l a n b n c m d , ^4 = C abcd n a m b n c m d . (38) 

When an electromagnetic field is present, we must include the complex tetrad components of 
the Maxwell field into the equations: 

00 = F ab l a m b , (39) 

01 = \F ab (I V + m a m b ) , 

02 = F ab n a m b , 

as well as the Ricci (or stress tensor) constructed from the three 0j, e.g., T ab l a l b — k(f>o(f> Q ,k = 
2Gc- 4 . 

Remark: We mention that much of the physical content and interpretations in the present 
work comes from the study of the lowest spherical harmonic coefficients in the leading terms of the 
far-field expansions of the Weyl and Maxwell tensors. 

The NP version of the vacuum (or Einstein-Maxwell) equations consists of three sets (or four 
sets) of nonlinear first-order coupled partial differential equations for the variables: the tetrad 
components, the spin coefficients, the Weyl tensor (and Maxwell field when present). Though there 



"7244, « = 7l31, 

"7243, P = 7134, 

-7242, O = 7133, 

"7241, T = 7132. 



(35) 
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is no hope that they can be solved in any general sense, many exact solutions have been found from 
them. Of far more importance, large classes of asymptotic solutions and perturbation solutions can 
be found. Our interest lies in the asymptotic behavior of the asymptotically-flat solutions. Though 
there are some subtle issues, integration in this class is not difficult [39, 45 . With no explanation 
of the integration process, except to mention that we use the Bondi coordinate and tetrad system 
of Equations ©, ©, and (fTTl) and asymptotic flatness, we simply give the final results. 

First, the radial behavior is described. The quantities with a zero superscript, e.g., a , ip®] • • • j 
are 'functions of integration', i.e., functions only of (ub,CiC)- 

• The Weyl tensor: 

ipo 

lj'2 
1p4 

• The Maxwell tensor: 

<t>i 

• The spin coefficients and metric variables: 

k = 7r = e = 0, r = a + /3, (42) 

p = p^-r- 1 -a°a°r- 3 + 0(r- 5 ), 

a = a°r- 2 + ((a°) 2 CT - ^ °/2) r- 4 + 0(r- 5 ), 

a = c/V- 1 +0(r~ 2 ), p = p°r- x + 0(r- 2 ), 

1 = 1° - ipUZr 2 )' 1 + 0(r~ 3 ), A = A°r~ 1 + 0(r~ 2 ), 

H = fjPr' 1 +0(r- 2 ), v = is Q + 0(r- 1 ), 





-5 _ 


VO{r 


- 6 ), 




-4 








-3 _ 




-% 




-2 


\-0(r 


-% 


<r~ 


-1 _ 


hO(K 





4^ 



0(r- 4 ) 
0(r- 3 ) 
0(r- 2 ) 



(40) 



(41) 



A 


= C or C, 










= e^- 1 - 




a a t 0A r- 3 -+ 


-0(r" 4 ), 


to 


= u/V" 1 - 


( ( 7 cJ + ^/2)r 


- 2 + 0(r- 3 


), 


X A 


= Wit 4 




+ 0(r- 4 ), 




u 


= f/°-(7° 


+ 7°)r-($ + 


^)(2 ? r 1 


+ 0(r- 2 ) 
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The functions of integration are determined, using coordinate conditions, as: 



s 


= -P £° C = 


f43) 


eOC 


= 0, £ = -P, 


( A A \ 

(44) 


D 
1 




(A Z\ 
(45) 


a 


P 2' 


(46) 


7° 


= u° = 0, 


(47) 




= -8ct°, 


(48) 


A° 


= cr , 


(49) 




= C7° = -l, 


(50) 


< 


-0 


(51) 




= OCT , 


(52) 




= g 2 CT°-g 2 CT + CT A -CT°A . 


(53) 



• The mass aspect, 

* = ^ 2 + g 2 CT + CT ^ , (54) 

satisfies the physically very important reality condition: 

* = (55) 

• Finally, from the asymptotic Bianchi identities, we obtain the dynamical (or evolution) rela- 
tions: 

iP 2 = -3^ + <rV2 + M$§, (56) 
= _g^° + 2ct°V3 + 2k<f)i$2, (57) 
4 = -g< + 3CTV 2 ° + 3^°, (58) 

(59) 
(60) 

k = 2Gc" 4 . (61) 

Remark: These last five equations, the first of which contains the beautiful Bondi energy- 
momentum loss theorem, play the fundamental role in the dynamics of our physical quantities. 

Remark: Using the mass aspect, \E', with Equations (fBTj) and (|5"2"1) . the first of the asymptotic 
Bianchi identities can be rewritten in the concise form, 

* = &W+k(f)^. (62) 

From these results, the characteristic initial problem can roughly be stated in the following 
manner. At ub = ubo we choose the initial values for (i/jq, ip®, ip®), i.e., functions only of (£, £). The 



A? = -30°, 

II = -3^ + ctV"; 
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characteristic data, the complex Bondi shear, <j°(ub,CiO> is then freely chosen. Since and ip® 
are functions of a , Equations ([49]) . (|5T|) and (|52|) and its derivatives, all the asymptotic variables 
can now be determined from Equations ((55)) - ([r?0j) . 

An important consequence of the SC formalism is that it allows simple proofs for many geo- 
metric theorems. Two important examples are the Goldberg-Sachs theorem [18] and the peeling 
theorem [ST]. The peeling theorem is essentially given by the asymptotic behavior of the Weyl 
tensor in Equation (|40[) (and Equation (|4Tj) ). The Goldberg-Sachs theorem essentially states that 
for an algebraically-special metric, the degenerate principle null vector field is the tangent field to 
a shear-free NGC. Both theorems are implicitly used later. 

One of the immediate physical interpretations arising from the asymptotically-flat solutions was 
Bondi's [5] identifications, at 3 + , of the interior spacetime four- momentum (energy /momentum). 
Given the mass aspect, Equation ([53]). 



and the spherical harmonic expansion 

* = *° + + >!'"> _,; , + . . . , (63) 

Bondi identified the interior mass and three-momentum with the I — and I — 1 harmonic contri- 
butions; 

Mb = (64) 

P< = -|L*<. (65) 

The evolution of these quantities, (the Bondi mass/momentum loss) is then determined from 
Equation ([r?2"]) . The details of this will be discussed in Section [SJ 

The same clear cut asymptotic physical identification for interior angular momentum is not as 
readily available. In vacuum linear theory, the angular momentum is often taken to be 

Jk = W Im { ^ k) - (66) 

However, in the nonlinear treatment, correction terms quadratic in er and its derivatives are of- 
ten included [55]. In the presence of a Maxwell field, this is again modified by the addition of 
electromagnetic multipole terms [35J [3] . 

In our case, where we consider only quadrupole gravitational radiation, the quadratic correction 
terms do in fact vanish and hence Equation ([66 j) . modified by the Maxwell terms, is correct as it is 
stated. 



2.5 The Bondi— Metzner— Sachs group 

The group of coordinate transformations at 3 + that preserves the Bondi coordinate conditions, the 
BMS group, is the same as the asymptotic symmetry group that arises from approximate solutions 
to Killing's equation as J + is approached. The BMS group has two parts: the homogeneous Lorentz 
group and the supertranslation group, which contains the Poincare translation sub-group. Their 
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importance to us lies in the fact that all the physical quantities arising from our identifications must 
transform appropriately under these transformations [49 , 29] . 

Specifically, the BMS group is given by the supertranslations, with a(£, Q an arbitrary regular 
diffcrcntiable function on S 2 : 

u B = mb + o;(C,C) (67) 

(C,C) = (CO 

and the Lorentz transformations, with (a, b, c, d) the complex parameters of SL(2, C), 

u B = Ku B , (68) 
1 + CC 



K = 



« + &)(3C + b) + « + d){c( + d) ' 



C = ad— be = 1. 

If a(C C) is expanded in spherical harmonics, 

a(c,c) = E am/1 w<^)< ( 69 ) 

the / = 0, 1 terms represent the Poincare translations, i.e., 

«(P) (C, = d a l a = y^d°Y ° - l -d l Y^. (70) 
Details about the representation theory, with applications, are given later. 

2.6 Algebraically-special metrics and the Goldberg— Sachs theorem 

Among the most studied vacuum spacetimes are those referred to as 'algebraically-special' space- 
times, i.e., vacuum spacetimes that possess two or more coinciding principle null vectors (PNVs). 
PNV fields [SD] (in general, four locally-independent fields exist) are defined by solutions, L a , to 
the algebraic equation 

L b L[ e C a ] bc [ d Lf]L c — 0, L a L a — 0. 
The Cartan-Petrov-Pirani-Penrose classification 52, 53, 50 describes the different degeneracies: 



Alg. General [1,1,1, 1] 

Type II [2,1,1] 

Type D or degenerate [2, 2] 

Type III [3, 1] 

Type IV or Null [4]. 

In NP language, if the tetrad vector l a is a principle null vector, i.e., L a = l a , then automatically, 
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For the algebraically-special metrics, the special cases are 

Type II Vo = ipi = 
Type III ipo = ipi = ip2 = 
Type IV ipa = ipi = ip2 = i'3 = 

rp n V'O = ^1 = ^3 = ^4 = 

±ype u with both / a and n a PNVs. 



An outstanding feature of the algebraically-special metrics is contained in the beautiful Goldberg- 
Sachs theorem |18j . 

Theorem. For a nonflat vacuum spacetime, if there is an NGC that is shear-free, i.e., there is a 
null vector field with (k = 0, a = 0), then the spacetime is algebraically special and, conversely, if a 
vacuum spacetime is algebraically special, there is an NGC with (k = 0, a = 0). 
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3 Shear-Free NGCs in Minkowski Space 

The structure and properties of asymptotically shear-free NGCs (our main topic) are best under- 
stood by first looking at the special case of congruences that are shear-free everywhere (except at 
their caustics). Though shear- free congruences are also found in algebraically-special spacetimes, 
in this section only the shear-free NGCs in Minkowski spacetime, M, are discussed [4] 



3.1 The flat-space good-cut equation and good-cut functions 

In Section [2J we saw that in the NP formalism, two of the complex spin coefficients, the optical 
parameters p and a of Equations (|27[) and (|28[) . play a particularly important role in their description 
of an NGC; namely, they carry the information of the divergence, twist and shear of the congruence. 

From Equations ([27| and (|28|) . the radial behavior of the optical parameters for general shear- free 
NGCs, in Minkowski space, is given by 

P = ^+ ^2 ■ CT = °> ( 71 ) 

where £ is the twist of the congruence. A more detailed and much deeper understanding of the 
shear-free congruences can be obtained by first looking at the explicit coordinate expression, Equa- 
tion (p3|. for all flat-space NGCs: 

x a = u B (l a + n a ) - Lm - Lrh a + (r* - r )f a , (72) 

where L(ub,00 is an arbitrary complex function of the parameters y w = (itB,C>C)i r 0j also an 
arbitrary function of (ubjCjC)i determines the origin of the affine parameter; r* can be chosen 
freely. Most frequently, to simplify the form of p, ro is chosen as 

r = -- (bL + BL + LL + Zl) . (73) 

At this point (|72l) describes an arbitrary NGC with £, () labeling the geodesies and r* the affine 
distance along the individual geodesies; later L(ub,(X) w iU be chosen so that the congruence is 
shear-free. 

The tetrad (l a , n a , m° , in ) is given by [27] 

x " '1 + CCC + CX-zC-i + CC), (74) 



2(1 + CO 

V2 
2(1 + CO 

V2 

2(1 + CO 



(i + CC,-(C + OX-<,i-CC), 
(o,i-cV*(i + C 2 ),2C). 



There are several important comments to be made about Equation (|72l) . The first is that there is a 
simple geometric meaning to the parameters (its, C) : they are the values of the Bondi coordinates 
of 3 + , where each geodesic of the congruence intersects 3 + . 

The second concerns the geometric meaning of L. At each point of 3 + , consider the past 
light cone and its sphere of null directions. Coordinatize that sphere (of null directions) with 
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stereographic coordinates. The function L(ub,CiC) * s the stereographic angle field on 3 + that 
describes the null direction of each geodesic intersecting 3 + at the point (ub,(,C)- values 
L = and L = oo represent, respectively, the direction along the Bondi l a and n a vectors. This 
stereographic angle field completely determines the NGC. 

The twist, E, of the congruence can be calculated in terms of L(ub,CjC) directly from Equa- 
tion (|72")l and the definition of the complex divergence, Equation (l2~i|) , leading to 



iS 



- |SL + LL-3L- Lt} . (75) 



We now demand that L be a regular function of its arguments, i.e., have no infinities, or, equivalently, 
that all members of the NGC come from the interior of the spacetime and not lie on 3 + itself. 

It has been shown [6] that the condition on the stereographic angle field L for the NGC to be 
shear-free is that 

(5L + LL = 0. (76) 

Our task is now to find the regular solutions of Equation (I76|) . The key to doing this is via the 
introduction of a new complex variable r and complex function [261 127] , 



r = T(u B X,(). (77) 

T is related to L by the CR equation (related to the existence of a CR structure on 3 + ; see 
Appendix [B]) : 

g (UB) r + LT = o. (78) 

Remark: The following 'gauge' freedom becomes useful later, r =>■ t* = F(t), with F analytic, 
leaving Equation (j78[) unchanged. In other words, 



r* = T* ( UB , C, C) = F (T(u B , C, 0) , (79) 

leads to 

g (ltB) T* = f3 (ttB) T, 
T* = F'f, 
d {uB) T*+LT* = 0. 

We assume, in the neighborhood of real 3 + , i.e., near the real ub and ( = C, that T(ub,CjC) 
is analytic in the three arguments (wb,C)0- The inversion of Equation (|77l) yields the complex 
analytic cut function 

u B = G(t,C,C)- (80) 

Though we are interested in real values for Mb, from Equation (|80p we see that for arbitrary r in 
general it would take complex values. Shortly, we will also address the important issue of what 
values of r are needed for real its- 

Returning to the issue of integrating the shear- free condition, Equation ([To]) , using Equa- 
tion (|77p . we note that the derivatives of T, 3( MB )T and T can be expressed in terms of the deriva- 
tives of G(t, £, C) by implicit differentiation. The ub derivative of T is obtained by taking the ub 
derivative of Equation (|5U|) : 

l = G'(r,C,Or^T=^-, (81) 
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while the <3(„ B )T derivative is found by applying S( MB ) to Equation (1801) . 

= G'(t,C,0^u b )T + B (t) G, (82) 



S(r)G 



G'(r,C,C) 

When Equations (|81[) and (|82|) are substituted into Equation ([78| , one finds that L is given implicitly 
in terms of the cut function by 

L(u B ,CX) = 3(r)G(r,C,C). (83) 
u B = G(t,C,C)^t = T(ub,C,C)- (84) 

Thus, we see that all information about the NGC can be obtained from the cut function 
G(t, CO- 

By further implicit differentiation of Equation ()83|) . i.e., 



s (ub) l(u b ,c,C) = 3[ r) G(T,c,c) + g (T) G'(T,(,().g (tlB) r, 
l(«b,C,0 = s M G'(r,c,C)-r, 

using Equation (|78|) , the shear-free condition ([76]) becomes 

9f r) G(r,C,C) = 0. (85) 

This equation will be referred to as the homogeneous Good- Cut Equation and its solutions as flat- 
space GCFs. In the next Sectional an inhomogeneous version, the Good- Cut Equation, will be 
found for asymptotically shear-free NGCs. Its solutions will also be referred to as GCFs. 

From the properties of the S 2 operator, the general regular solution to Equation (|85[) is easily 
found: G must contain only I = and I = 1 spherical harmonic contributions; thus, any regular 
solution will be dependent on four arbitrary complex parameters, z a . If these parameters are 
functions of r, i.e., z a = £ a (r), then we can express any regular solution G in terms of the complex 
world line £°(t) [2S[27]: 

u B = G(t, C, C) = e(r)la(CX) = - (86) 

The angle field L(ub, C; C) then has the form 

L(u B ,C,0 = S(r)G(r,C,C)=r(r)m (C,C). (87) 
ub = C(r)ia(CX)- 



Thus, we have our first major result: every regular shear-free NGC in Minkowski space is 
generated by the arbitrary choice of a complex world line in what turns out to be complex Minkowski 
space. See Equation ([70)) for the connection between the I = (0, 1) harmonics in Equation (l86l) and 
the Poincare translations. We see in the next Section [4] how this result generalizes to regular 
asymptotically shear-free NGCs. 
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Remark: We point out that this construction of regular shear-free NGCs in Minkowski space 
is a special example of the Kerr theorem (c.f., [51]). Writing Equations (|87l) and ([88]) as 

a + b( + 6C + cCC 

1 + CC ' 

(b + cQ-ga + bC) 

i + CC 

where the (a(r), b(r), c(r), rf(r)) are simple combinations of the £°(r), we then find that 

L + uC = b + c(, 
u — LQ = a + b(. 

Noting that the right-hand side of both equations are functions only of r and we can eliminate 
the t from the two equations, thereby constructing a function of three variables of the form 

F{L + u(,u- L(,() = 0. 

This is a special case of the general solution to Equation (|75|) . which is the Kerr theorem. 

In addition to the construction of the angle field, L(itB,C>C)> from the GCF, another quantity 
of great value in applications, obtained from the GCF, is the local change in u B as t changes, i.e., 

V{T,C,l)=d T G = G'. (89) 
3.2 Real cuts from the complex good cuts, I 

Though our discussion of shear-free NGCs has relied, in an essential manner, on the use of the 
complcxification of 3 + and the complex world lines in complex Minkowski space, it is the real 
structures that are of main interest to us. We want to find the intersection of the complex GCF 
with real 3 + , i.e., what are the real points and real cuts of ub = G(r, £, £), (C — C)> an d what are 
the values of r that yield real ub- 

To construct an associated family of real cuts from a GCF, we begin with 

u B = G(r,C,C) = ^e°(r) - if (r)r°(C,C) (90) 

and write 

r = s + iX (91) 
with s and A real. The cut function can then be rewritten 

u B = G(r,C,C) =_G(s + iX, CO _ (92) 
= Gr(s, A, C, C) + iGi(s, A, C)C)i 

with real Gr(s 7 A, £, () and G/(s, A, C, C)- The Gr(s, A, C, C) an d G/(s, A, (, £) are easily calculated 
from G(t, C, C) by 

G R (s,X,(,0 = ^{G(s + i\,(,0 + G(s + i\, (,()}, (93) 
Gj(s,XXX) = ^{G(s + iX,(,0-G(s + iX, (,()}■ (94) 



U B = 

L = 
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By setting 
and solving for 

we obtain the associated real slicing, 



G J («,A,C,C) = 
A = A(s,C,C) 



u { B R) =G R (s,A(s, c,0» CO- 
Thus, the values of r that yield real values of ub are given by 

r = s + iA(s,CO- 
As an example, using Equation (|90[) . we find to first order in A 



,(«) 



.1 



^(C,0 



^(CO, 



G fl ( S ,A,C,0 



A(C,0 = -i 



^?(s) + 0(s)r°(c,0 



n°i(C,0, 



[V2&(*)' 01' 

An Important Remark: We saw earlier that the shear-free angle field was given by 

L(u B ,C.O = 3(r)G(r,C,0, 

UB = G(t,C,0^t = T( Ub ,C,0, 



(95) 
(96) 

(97) 
(98) 

(99) 
(100) 
(101) 



(102) 
(103) 



where real values of ub should be used. If the real cuts, ub — Gr (s, A(s, Q, 0, C) > were used 
instead to calculate L(ub,00> the results would be wrong. The restriction of r to yield real ub, 
does not commute with the application of the 3 operator, i.e., 

L(KB,C,0^9G fi . 

The 3 differentiation must be done first, holding r constant, before the reality of ub is used. In other 
words, though we are interested in real 3 + , it is essential that we consider its (local) complexification. 



3.3 Approximations 

Due to the difficulties involved in the intrinsic nonlinearities and the virtual impossibility of exactly 
inverting arbitrary analytic functions, it often becomes necessary to resort to approximations. The 
basic approximation will be to consider the complex world line ^°(t) as being close to the straight 
line, £o( T ) = tSq] deviations from this will be considered as first order. We retain terms up to 
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second order, i.e., quadratic terms. Another frequently used approximation is to terminate spherical 
harmonic expansions after the I = 2 terms. 

It is worthwhile to discuss some of the issues related to these approximations. One important 
issue is how to use the gauge freedom, Equation (|79|). r => r* = <£>(t), to simplify the 'velocity 
vector', 



«»(r) = e'(r) = (104) 

A Notational issue: Given a complex analytic function (or vector) of the complex variable r, 
say G(r), then G(t) can be decomposed uniquely into two parts, 

G(t) = <5 R (T) + i<5i(T), 

where all the coefficients in the Taylor series for <8r(t) and <5i(r) are real. With but a slight 
extension of conventional notation we refer to them as real analytic functions. 
With this notation, we also write 

By using the reparametrization of the world line, via r* = $(t), with the decomposition 

$(t) = $ fl (r) +i&i(r), 

the 'velocity' transforms as 

v* a (r*) =w a (r)[$(r)']- 1 . 

One can easily check that by the appropriate choice of $j(t) we can make v R a (T*) and v} a (r*) 
orthogonal, i.e., 

VabVRvf = 0, (105) 

and by the choice of $^(r), the v*^(t*) can be normalized to one, 

Va b v* R a v* R b = 1. (106) 

The remaining freedom in the choice of $(r) is simply an additive complex constant, which is 
used shortly for further simplification. 

We now write v r (t) — (v r (t), v r (t)^ , which, with the slow motion approximation, yields, from 
the normalization, 

<(r) = ^l + {v R Y « 1 + \vf + ■ ■ ■ (107) 
From the orthogonality, we have 

«iVWh(tK(t), 

i.e., Vj(t) is second order. Since U°(r) = C?( r )' i s second order, ^j(r) is a constant plus a second- 
order term. Using the remaining complex constant freedom in $(t), the constant can be set to 
zero: 

£°( T ) = second order. (108) 
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Finally, from the reality condition on the ub, Equations (|94p . (|§T|) and (IM1) yield, with t = s + «A 
and A treated as small, 



4f } = &(s%+vtts)i a -l^, (109) 
A = A( S ,C,C) = -^ £ ^, (HO) 

i-^'w; 

Within this slow motion approximation scheme, we have from Equations (|109l) and (IllOp . 

vg? = V2u^ =s-^Us)Y 1 i +2vUs)i a e i (s)i b , (111) 

A « (H2) 
or, to first order, which is all that is needed, 



We then have, to linear order, 



r = s + i^Q(s)Y°, (113) 



2 



3.4 Asymptotically-vanishing Maxwell fields 
3.4.1 A prelude 

The basic starting idea in this work is simple. It is in the generalizations and implementations 
where difficulties arise. 

Starting in Minkowski space in a fixed given Lorentzian frame with spatial origin, the electric 
dipole moment ~t) e is calculated from an integral over the (localized) charge distribution. If there 
is a shift, % in the origin, the dipole transforms as 

Di = l5 E -qt. (114) 

If i3 e is time dependent, we obtain the center-of-charge world line by taking D* E = 0, i.e., from 
it = ~D eQ 1 - It is this idea that we want to generalize and extend to gravitational fields. 
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The first generalization is formal and somewhat artificial: shortly it will become quite natural. 
We introduce, in addition to the electric dipole moment, the magnetic dipole moment ~~t) m (also 
obtained by an integral over the current distribution) and write 

~i>c = D E +W M - 

By allowing the displacement it to take complex values, i^c, Equation (| 1 14[) . can be generalized 
to _^ 

D* c =t c - q%, (115) 

so that the complex center-of-charge is given by 15^ — or 

~R\ c = ~3 c q-\ (116) 

We emphasize that this is done in a fixed Lorentz frame and only the origin is moved. In different 
Lorentz frames there will be different complex centers of charge. 

Later, directly from the general asymptotic Maxwell field itself (satisfying the Maxwell equa- 
tions), we define the asymptotic complex dipole moment and give its transformation law, including 
transformations between Lorentz frames. This yields a unique complex center of charge independent 
of the Lorentz frame. 



3.4.2 Asymptotically-vanishing Maxwell fields: General properties 

Our first set of applications of shear-free NGCs comes from Maxwell theory in Minkowski space. We 
review the general theory of the behavior of asymptotically-flat or vanishing Maxwell fields assuming 
throughout that there is a nonvanishing total charge, q. As stated in Section the Maxwell field 
is described in terms of its complex tetrad components, (4>o,4>i,<j>2)- In a Bondi coordinate/tetrad 
system the asymptotic integration is relatively simple [35l [25] resulting in the radial behavior (the 
peeling theorem): 



1 0(r- 4 ), (117) 



r 

where the leading coefficients of r, (0°, 0°, (j)®) satisfy the evolution equations: 

$ + 8$ - 0, (118) 
$ + 3$ = 0. (119) 

The formal integration procedure is to take </>2 as an arbitrary function of (mb,CjC) (the free 
broadcasting data), then integrate the second, for <ftl, with a time- independent spin- weight s = 
function of integration and finally integrate the first, for <pQ. Using a slight modification of this, 
namely from the spherical harmonic expansion, we obtain, 

= ^i + <i y 2« + ---. (120) 

4>1 = q + 4> uY 1 i + 4> HjY° j + ..., (121) 
4>l = ^i^+^yJW + •••> ( 122 ) 
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with the harmonic coefficients related to each other by the evolution equations: 

4 = 2qr 1 i (u let )Y 1 \ + Q i i%\ j + ..., (123) 

$ = -2gr, i »(i heb )Y^-^' u Y^ + ... 
The physical meaning of the coefficients are 

q = total electric charge, (124) 
qrf — D l c — complex (electric & magnetic) dipole moment = D l E + iD l M , 
= complex (electric & magnetic) quadrupole moment, 

etc. For later use, the complex dipole is written as Dp(it rc t) = qrf(u lc t). Note that the D z c is 
defined relative to a given Bondi system. This is the analogue of a given origin for the calculations 
of the dipole moments of Equation (|114[) . 

Later in this section it will be shown that we can find a unique complex world line, ^"(t) = 
(the world line associated with a shear-free NGC), that is closely related to the rj l (u let ). 
From this complex world line we can define the intrinsic complex dipole moment, D l xc — q£?(s). 

However, we first discuss a particular Maxwell field, F ab , where one of its eigenvectors is a 
tangent field to a shear-free NGC. This solution, referred to as the complex Lienard-Wiechert field 
is the direct generalization of the ordinary Lienard-Wiechert field. Though it is a real solution in 
Minkowski space, it can be thought of as arising from a complex world line in complex Minkowski 
space. 

3.4.3 A coordinate and tetrad system attached to a shear-free NGC 

The parametric form of the general NGC was given earlier by Equation (j72|) . 

x a = u B (l a + n a ) - Dm - Lm a + (r* - r )l a . (125) 

The parameters (ub, Cj C) labeled the individual members of the congruence while r* was the affinc 
parameter along the geodesies. An alternative interpretation of the same equation is to consider it 
as the coordinate transformation between the coordinates, x a (or the Bondi coordinates) and the 
geodesic coordinates (ub, r*, C)- Note that these coordinates are not Bondi coordinates, though, 
in the limit, at 3 + , they are. The associated (geodesic) tetrad is given as a function of these geodesic 
coordinates, but with Minkowskian components by Equations (f74|) . We restrict ourselves to the 
special case of the coordinates and tetrad associated with the L from a shear- free NGC. Though 
we are dealing with a real shear-free twisting congruence, the congruence, as we saw, is generated 
by a complex analytic world line in the complexified Minkowski space, z a = £ (t). The complex 
parameter, r, must in the end be chosen so that the 'its' of Equation (|90| is real. The Minkowski 
metric and the spin coefficients associated with this geodesic system can be calculated [27] in the 
(u-B,r* , C) C) frame. Unfortunately, it must be stated paramctrically, since the r explicitly appears 
via the £ a (r) and can not be directly eliminated. (An alternate choice of these geodesic coordinates 
is to use the r instead of the ub- Unfortunately, this leads to an analytic flat metric on the 
complexified Minkowski space, where the real spacetime is hard to find.) 
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The use and insight given by this coordinate/tetrad system is illustrated by its application to a 
special class of Maxwell fields. We consider, as mentioned earlier, the Maxwell field where one of its 
principle null vectors, l* a , (an eigenvector of the Maxwell tensor, F a t,l* a — \l%), is a tangent vector 
of a shear- free NGC. Thus, it depends on the choice of the complex world line and is therefore 
referred to as the complex Lienard-Wiechert held. (If the world line was real it would lead to the 
ordinary Lienard-Wiechert held.) We emphasize that though the source can formally be thought 
of as a charge moving on the complex world line, the Maxwell held is a real held on real Minkowski 
space. It will have a real (distributional) source at the caustics of the congruence. Physically, its 
behavior is very similar to real Lienard-Wiechert helds, the essential difference is that the electric 
dipole is now replaced by the combined electric and magnetic dipoles. The imaginary part of the 
world line determines the magnetic dipole moment. 

In the spin-coefficient version of the Maxwell equations, using the geodesic tetrad, the choice of 
l* a as the principle null vector 'congruence' is just the statement that 

<f>l = F ab l* a m* b = 0. 

This allows a very simple exact integration of the remaining Maxwell components [35] . 



3.4.4 Complex Lienard Wiechert Maxwell field 

The complex Lienard-Wiechert fields (which we again emphasize are real Maxwell fields) are for- 
mally given by the (geodesic) tetrad components of the Maxwell tensor in the null geodesic coor- 
dinate system (ub , r* , £, Q , Equation (|125[) . As the detailed calculations are long [35] and take us 
too far afield, we only give an outline here. The integration of the radial Maxwell equations leads 
to the asymptotic behavior, 



0, (126) 
P 2 K, (127) 
p^ + O(p 2 ), (128) 



with 



P = -(r*+iE)-\ (129) 
2iH = (5L + LL — 9Z — LL. 

The order expression is known in terms of (</>*°, </>2°)- The function Hub, £, £) is given by 

L{ubXX) = 3(t)G(t,(,C), 

u B = G(r,(,(>r(r)f„((,(), 

with £ a (r) an arbitrary complex world line that determines the shear- free congruence whose tangent 
vectors are the Maxwell held eigenvectors. 

Remark: In this case of the complex Lienard-Wiechert Maxwell field, the £ a determines the 
intrinsic center-of-charge world line, rather than the relative center-of-charge line. 

The remaining unknowns, 01°, </>2°, are determined by the last of the Maxwell equations, 

B# +Z&°+2£# = 0, (130) 
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which have been obtained from Equations (|5§)) and (pD|) via the null rotation between the Bondi 
and geodesic tetrads and the associated Maxwell field transformation, namely, 

l a -> l* a = l a - -m a - -fh a + 0(r*- 2 ), (131) 
r r 

m a -> m* a = m a — -n a , (132) 
r 

n a -> n* a = n a , (133) 



with 



61° - " 

A*0 J,0 

»2 = 02 



O = 0°-2L0? + L 2 0°, (134) 
? - (135) 

(136) 



These remaining equations depend only on L(ub, Cj C)j which, in turn, is determined by £ a (r). In 
other words, the solution is driven by the complex line, £°(t). As they now stand, Equations (|130[) 
appear to be difficult to solve, partially due to the implicit description of the L(ub, C C)- 

Actually they are easily solved when the independent variables are changed, via Equation 
from (mb, C, C) to the complex (r, £, £). They become, after a bit of work, 



3(r)(^ 2 0?) - 0, (137) 
3( T )M) = 0?', (138) 
V = r'(r)L(C,C), (139) 

with the solution 

4>T = §^ 2 , (140) 

(7 being the Coulomb charge. 

Though we now have the exact solution, unfortunately it is now in complex coordinates where 
virtually every term depends on the complex variable r, via £ a ( r )- This is a severe impediment to 
a full description and understanding of the solution in the real Minkowski space. 

In order to understand its asymptotic behavior and physical content, one must transform it, 
via Equations (|131|) - (|136j) . back to a Bondi coordinate/tetrad system. This can only be done by 
approximations. After a lengthy calculation [35], we find the Bondi peeling behavior 

00 = r- 3 ^ + 0(r- 4 ), (141) 

01 = r- 2 ^ + 0(r- 3 ), 
4> 2 = r- J ^ + 0(r- 2 ), 
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g = q[LV-' 1 + ^L 2 V-% T) [V- l V'}), (142, 
' i = ^(1 + LVB (T) [V- 1 V']), (143) 



with 

1 
2 

i.o _ 1 ft i n/a. . n/-ii/' 
2U 2 

$ - -f^-^M^-V], (144) 

v = r'uoc). (145) 

Next, treating the world line, as discussed earlier, as a small deviation from the straight line, 
e(r) = r*g, i.e., by 

CM - (r + ^°(r),r(r)), 
£*( r ) <C 1, <5£° (r) = second order. 
The GCF and its inverse (see Section [5]) are given, to first order, by 

u rct = V2u B = V2G = T-^f(r)yO(C,C), (146) 

T = "ret + ^fKet)^(C,0- (147) 

Again to first order, Equations (|T42|) . (fT43|) and (|T44|) yield 

0° = 2qC(u 1 . ct )Y 1 \, (148) 

= q + V2qe'(u let )Y°, 
& = -2qe"(u rct )Y u \ 

the known electromagnetic dipole field, with a Coulomb charge, q. One then has the physical 
interpretation of £ a (u ro t) : q^ a (u Te t) is the complex dipole moment; (the electric plus times 
magnetic dipole) and £ a (u rct ) is the complex center of charge, the real part being the ordinary 
center of charge, while the imaginary part is the 'imaginary' magnetic center of charge. This simple 
relationship between the Bondi form of the complex dipole moment, g^*(ti ret ), and the intrinsic 
complex center of charge, ^"(t), is true only at linear order. The second-order relationship is given 
later. 

Reversing the issue, if we had instead started with an exact complex Lienard-Wiechert field but 
now given in a Bondi coordinate/tetrad system and performed on it the transformations, Equa- 
tions (fT4")) and (|134p to the geodesic system, it would have resulted in 

05 = o. 

This example was intended to show how physical meaning could be attached to the complex 
world line associated with a shear- free NGC. In this case and later in the case of asymptotically-flat 
spacetimes, when the GCF is singled out by either the Maxwell field or the gravitational field, it 
will be referred to it as a UCF. For either of the two flat-space asymptotically- vanishing 

Maxwell field (with nonvanishing total charge) and for a vacuum asymptotically-flat spacetime, 
there will be a unique UCF. In the case of the Einstein-Maxwell fields there, in general, will be two 
UCFs, one for each field. 
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3.4.5 Asymptotically vanishing Maxwell fields & shear-free NGCs 

We return now to the general asymptotically- vanishing Maxwell field, Equations (|117|) and (|120[) . 
and its transformation behavior under the null rotation around n a , 

l a -> l* a = l a - -m a - ~m° + 0(r~ 2 ), (149) 
r r 

m a -> m* a =m a --n a , 
r 

rt — > n = n , 

with L(ub,OC) = ■? a ( r ) r "'aj being one of our shear-free angle fields defined by a world line, z a = 
£ a (r). The leading components of the Maxwell fields transform as 



L*0 lO 



bl-2L(j^ + L 2 <jP 2 , (150) 



- (151) 

(152) 



6*° = 



The 'picture' to adopt is that the new (j)*s are now given in a tetrad defined by the complex 
light cone (the generalized light cone) with origin on the complex world line. (This is obviously 
formal and physically nonsense, but mathematically quite sound, as the shear-free congruence can 
be thought of as having its origin on the complex line, £ q (t).) From the physical identifications 
of charge, dipole moments, etc., of Equation (|123[) . we can obtain the transformation law of these 
physical quantities. In particular, the I — 1 harmonic of </>q, or, equivalently, the complex dipole, 
transforms as 

C = C-2(i0 1 )| i . + (i 2 ^)k (153) 

where the notation W\i means extract only the I — 1 harmonic from a Clebsch-Gordon expansion 
of W . A subtlety and difficulty of this extraction process is explained/clarified below. 

If by some accident the Maxwell field was a complex Lienard-Wiechert field, a world line £°(t) 
could be chosen so that 0q° = 0- However, though this cannot be done in general, the / = 1 
harmonics of can be made to vanish by the appropriate choice of the £ a (r). This is the means 
by which a unique world line is chosen. 

A Clarifying Comment: 

An important observation, obvious but easily overlooked, concerning the spherical harmonic ex- 
pansions, is that, in a certain sense, they lack uniqueness. As this issue is significant, its clarification 
is important. 

Assume that we have a particular spin-s function on 3 + , say, ?7( s )(ltB, C> C)i given in a specific 
Bondi coordinate system, (ttB,C>C)) that has a harmonic expansion given, for constant «b, by 

r ?(s) ( U B,C,C) = i(j | J^Hu^Y^ 

If exactly the same function was given on different cuts or slices, say, 

ub = G(s,(,0, (154) 

with 

»?W («» CC) = »?(.) (G( s ,c,C),C,C), 
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the harmonic expansion at constant s would be different. The new coefficients are extracted by the 
two-sphere integral taken at constant s: 

v *Um...) (s)= I ^(sX^U.^S. (155) 
J s 2 

It is in this rather obvious sense that the expansions are not unique. From Equation ()97[) . using 
the real values of the complex cuts, i.e., 

u B = G R (s,A(s,t,0,C,Z), (156) 

(in place of u-q = G{s, C)C))i the transformation, Equation (|153|) . and harmonic extraction yields 
(f>Q*(s). This is proportional to the transformed complex dipole moment. This transformation, a 
functional of the form, 

08?(*)= FiC^US, (157) 
is decidedly nontensorial: in fact it is very nonlocal and nonlinear. 

Though it is clear that extracting (%* (s) with this relationship is available in principle, in practice 
it is impossible to do it exactly and all examples are done with approximations: essentially using 
slow motion for the complex world line. 

3.4.6 The complex center of charge 

The complex center of charge is defined by the vanishing of the complex dipole moment (s); in 
other words, 

ri($US,$,£°) = o (158) 

determines the (up to this point) arbitrary complex world line, £ q (t). In practice we do this only 
up to second order with the use of only the (I — 0, 1,2) harmonics. The approximation we are 
using is to consider the charge q as zeroth order and the dipole moments and the spatial part of 
the complex world line as first order. 
From Equation (I153j) . 

C =r < «0§ < -2L4\ i = o (159) 

with the identifications, Equation (|123p . for q and Dc, we have to first order (with s ~ u re % ~ r), 

£>cKet) = qrf{u TCt ) = qC(u let ). (160) 

This is exactly the same result as we obtained earlier in Equation (II , via the charge and current 
distributions in a fixed Lorentz frame. 

Carrying this calculation |35] to second order, we find the second-order relationship between the 
intrinsic complex dipole and the intrinsic dipole, 

Dl c = qC(s), D i c = q V i (s), (161) 

e( S ) = V k (s) + ^V i^)v k/ (s)~ l -e ljk ^(s)^( S ). (162) 

In Section [5l these ideas are applied to GR, with the complex electric and magnetic dipoles 
being replaced by the complex combination of the mass dipole and the angular momentum. 
The GCF, 

«B=e°(r)ia(C,C). 

with this uniquely determined world line is referred to as the Maxwell UCF. 
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4 The Good-Cut Equation and %-Space 

In Section [3l we discussed NGCs in Minkowski spacetime that were shear-free. In this section 
we consider asymptotically shear-free NGCs in asymptotically-flat spacetimes. That is to say, we 
consider NGCs that have nonvanishing shear in the interior of the spacetime but where, as null 
infinity is approached, the shear vanishes. In this case, whereas fully shear-free NGCs almost 
never occur, asymptotically shear-free congruences always exist. The case of algebraically-special 
spacetimes is the exception; they do allow one or two shear-free congruences. 

We begin by reviewing the shear-free condition and follow with its generalization to the asymp- 
totic shear-free case. From this we derive the generalization of the homogeneous good-cut equation 
to the inhomogeneous good-cut equation. Almost all the properties of the shear-free and asymp- 
totically shear-free NGCs come from the study of these equations. We will see that virtually all 
the attributes of shear-free congruences are shared by the asymptotically shear-free congruences. 
It is from the use of these shared attributes that we will be able to extract physical identifications 
and information (e.g., complex center of mass/charge, Bondi mass, linear and angular momentum, 
equations of motion, etc.) from the asymptotic gravitational fields. 

Though again the use of the formal complexification of 3 + , 1.6.; 5 IS essential for our analysis, 
it is the extraction of the real structures that is important. 

4.1 Asymptotically shear-free NGCs and the good-cut equation 

We saw in Section |3] that shear- free NGCs in Minkowski space could be constructed by looking at 
their properties near 3 + , in one of two equivalent ways. The first was via the stereographic angle 
field, L(ub,CiC)i which gives the directions the null rays make at their intersection with 3 + . The 
condition for the congruence to be shear-free was that L must satisfy 

d (uB) L + LL = 0. (163) 

We required solutions that were all nonsingular (regular) on the (£, £) sphere. (This equation has 
in the past most often been solved via twistor methods [19].) 

The second was via the complex cut function, ub = G(t, £, 0, that satisfied 

df T) G = 0. (164) 

The regular solutions were easily given by 

u B = G(T,c,o = r(T)wc,c) (165) 

with inverse function, 

t = r(u B ,C>0- 

They determined the L(ub,CiC) that satisfies Equation (|163p by the parametric relations 

L(u B XX) = S (T )G(r,C,C), (166) 
ub = r(rK(C,C), 

or by 

L{u B , C, C) = 3(t)G(t, C, C)I t =t(«b,C,C)' 
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where £ a (r) was an arbitrary complex world line in complex Minkowski space. 

It is this pair of equations, (|163[) and (|164[) . that will now be generalized to asymptotically-flat 
spacetimes. 

In Section [21 we saw that the asymptotic shear of the (null geodesic) tangent vector fields, l a , 
of the out-going Bondi null surfaces was given by the free data (the Bondi shear) o-°(mb, C)- If, 
near 3 + , a second NGC, with tangent vector l* a , is chosen and then described by the null rotation 
from l a to l* a around n a by 

l* a = l a + bm a + bm a + bbn a , (167) 
n* a = m 
n* a = n a , 

b = -L/r + 0(r- 2 ), 

with L(ub, Cj C) an arbitrary stereographic angle field, then the asymptotic Weyl components trans- 
form as 

L 3 r 3 °+L^l (168) 
bl (169) 

(170) 
(171) 
(172) 

and the (new) asymptotic shear of the null vector field l* a is given by [6l [26] 

a * =a° -d {uB) L-LL. (173) 

By requiring that the new congruence be asymptotically shear-free, i.e., a * — 0, we obtain the 
generalization of Equation (|163[) for the determination of L(ub, Cj C)j namely, 

Q {uB) L + LL = a°(u B ,(X)- (174) 

To solve this equation we again complexify 3 + to by freeing C to C an d allowing ub to take 
complex values close to the real. 

Again we introduce the complex potential r = T(ub, Cj C) that is related to L by 

S (ub) T + LT = 0, (175) 

with its inversion, 

ub = G(t,C,0- (176) 
Equation (|174|) becomes, after the change in the independent variable, ub =>■ r = T(wbjCjC)j an< i 



r ° 


= <- 


4LVi*° 


+ 6L 2 ?P* 




= 


3L^2 4 








2L^ 4 






= 






02 









implicit differentiation (see Section I37LI for the identical details), 

g^ T) G = <7°(G,C,C)- (177) 
This, the inhomogeneous good-cut equation, is the generalization of Equation (|164[) . 
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In Section 14.21 we will discuss how to construct solutions of Equation (|177p of the form, ub — 
G(t, C)C)j however, assuming we have such a solution, it determines the angle field L(ub, (, C) by 
the parametric relations 

L(« B> C,0 = 3<t)G, (178) 
u B = G(r,C,(). 



We return to the properties of these solutions in Section 14.21 

4.2 H-space and the good-cut equation 

Equation (|177|) . written in earlier literature as 

<5 2 Z = a°(ZXX), (179) 

is a well-known and well-studied partial differential equation, often referred to as the "good-cut equa- 
tion" [T9l [20] . For sufficiently regular <j (ub,(iO (which is assumed here) it has been proven [20] 
that the solutions are determined by points in a complex four-dimensional space, z a , referred to as 
H-space, i.e., solutions are given as 

u B = Z(z a ,(,0- (180) 

Later in this section, by choosing an arbitrary complex analytic world line in "H-space, z a = 
£ q (t), we describe how to construct the shear-free angle field, L(ub, C, 0- First, however, we discuss 
properties and the origin of Equation (|180[) . 

Roughly or intuitively one can see how the four complex parameters enter the solution from the 
following argument. By writing Equation (I179p as the integral equation 



Z = z Q Z a (C, + f <r°(Z, r], v)K+_ 2 ( Vl fj, C, OdS v (181) 

with 

i (i + C^ 2 (^-C) 



4tt (l + COCl+wiW-C)' 



'(i + w) 2 ' 



where z a l a ((, £) is the kernel of the S 2 operator (the solution to the homogeneous good-cut equation) 
and Kq _ 2 (C, Ci Vi V) 1S the Green's function for the <3 2 operator [23]. By iterating this equation, 
with the kernel being the zeroth iterate, i.e., 



Z„(C,C) = z a UCX)+ K+_ 2 (C,C,v,v)a(Z n ^,r),fj)dS v , (182) 

Js 2 

Z (C,0 - W„(C,C), (183) 

one easily sees how the four z a enter the solution. Basically, the z a come from the solution to the 
homogeneous equation. 
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It should be noted again that the z a l a ((, Q is composed of the I = (0, 1) harmonics, 

zX(CC) = ^=z° - ^(CO- (184) 

Furthermore, the integral term does not contribute to these lowest harmonics. This means that 
solutions can be written 

u B = Z(z a ,(, = z a UC, + Z t > 2 (z a , c, C), (185) 

with Zi>2 containing spherical harmonics I = 2 and higher. 

We note that using this form of the solution implies that we have set stringent coordinate 
conditions on the "H-space by requiring that the first four spherical harmonic coefficients be the 
four "H-space coordinates. Arbitrary coordinates would just mean that these four coefficients were 
arbitrary functions of other coordinates. How these special coordinates change under the BMS 
group is discussed later. 

Remark: It is of considerable interest that on 'H-space there is a natural quadratic complex 
metric - demonstrated in Appendix |D] - that is given by the surprising relationship |34[ 120] 

dsf H) = 9 { n )ab dz«dz^(±-J s ^) \ (186) 
dZ = V a Zdz a , (187) 
dS = ^^L. (188) 

(i + CO 2 

Remarkably this turns out to be a Ricci-fiat metric with a nonvanishing anti-self-dual Weyl tensor, 
i.e., it is intrinsically a complex vacuum metric. For vanishing Bondi shear, H-space reduces to 
complex Minkowski space (i.e., g(H)ab\<j°=o = Vab)- 

4.2.1 Solutions to the shear-free equation 

Returning to the issue of the solutions to the shear-free condition, i.e., Equation (|174[l . L[u B , C> C)j 
we see that they are easily constructed from the solutions to the good-cut equation, ub = Z(z a , £, (). 
By choosing an arbitrary complex world line in the "H-space, i.e., 

z a =a a (T), (189) 

we write the GCF as 

u B =G(r,C,0 = Z(e(r),C,0, (190) 

or, from Equation (|185|) . 

u B = G(t, C, C) = ^C°(r) - if (r)yO(C, Q + C j (r)Y 2 %(C, £) + •••• (191) 

This leads immediately, via Equations (|178[) and (| 1 79[) , to the parametric description of the shear- 
free stereographic angle field L(u B , (, C)j as wei l as the Bondi shear <j°(u B , Ci C) : 

y|C ( T ) — 2 ^ 1 ' ' ' i - ' ^ • ^ ' r ^ 1 ' " 2, i 

L(ub,(X) = f(T)Y 1 1 2 (C,C)-6^(r)r 2 ^.(C,C) + ---, (193) 
*°(u B ,U) = 24C j (r)Y 2 % + .... (194) 



u B = _e°(r)--e(r)F 1 l (C,C)+r J (r)F 2 ° I? (C,C) + ---, (192) 
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We denote the inverse to Equation (|19ip by 

T = r(u B ,C,C)- (195) 

The asymptotic twist of the asymptotically shear-free NGC is exactly as in the flat-space case, 

is = -|gr+ii-g£-ii|. (196) 

As in the flat-space case, the derived quantity 

V[r,C,l) = d T G = G' (197) 

plays a large role in applications. (In the case of the Robinson- Trautman metrics i55, 28] V is the 
basic variable for the construction of the metric.) 

Using the gauge freedom, r => r* = $(t)), in a slightly different way than in the Minkowski- 
space case, we impose the simple condition 

C° = t. (198) 

A Brief Summary: The description and analysis of the asymptotically shear-free NGCs in 
asymptotically-flat spacetimes is remarkably similar to that of the flat-space regular shear-free 
NGCs. We have seen that all regular shear-free NGCs in Minkowski space and asymptotically-flat 
spaces are generated by solutions to the good-cut equation, with each solution determined by the 
choice of an arbitrary complex analytic world line in complex Minkowski space or H-space. The 
basic governing variables are the complex GCF, u-q = G(t, C>C); an d the stereographic angle field 
on , L(ub, Ci C)> restricted to real 3 + . In every sense, the flat-space case can be considered as a 
special case of the asymptotically-flat case. 

In Sections [5] and [6l we will show that in every asymptotically flat spacetime a special complex- 
world line (along with its associated NGC and GCF) can be singled out using physical considera- 
tions. This special GCF is referred to as the (gravitational) UCF, and is denoted by 

u b =X(t,(,Q. (199) 
4.3 Real cuts from the complex good cuts, II 

The construction of real structures from the complex structures, i.e., finding the complex values of 
r that yield real values of ub and the associated real cuts, is virtually identical to the flat-space 
construction of Section [3] The only difference is that we start with the GCF 

u B = G(r, C,C) = r(r)l(C, + G ; > 2 (r, £ C) (200) 

rather than the flat-space 

u B = G(T,CX)=e(r)i a (C,0- 
Using t = s + i\, Equation (|200[) is written 

«4 B) =G R (s,\CO+iGi(s ) X > C,0- (201) 
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The reality of ub, i.e. 

again leads to 

and the real slicing, 



G J ( S) A,C,C) = 0, 
A = A( S ,C,C) 



,(R) _ 



= G R (s,A( s ,c,c),C,C) 



Using Equation (|191[) and expanding to first order in A, we have the expressions: 



V2 



\&(s) - \e i {s) , X--[CH{s)-Q{s) , X]Y^ i (C > 



(202) 
(203) 
(204) 

(205) 



+[^)-e/( S yx]Y 2 %(cx) 



+i 



G R (s,A(sXX)X,C) 

+[ei{s)-ei{syx]Y^{cro, 



(206) 



A = A( S ,C,C) = - 



(207) 



Since in all applications A(s, £, C) is multiplied by a first-order quantity, only the first-order expres- 
sion for A(s, £, £) is needed. Using the slow motion assumption, we have, 



a = ^K*)n°i(c,c) - v^>)f4(c,o. 



(208) 



Substituting this into Equation ()206[) , via the slow motion assumption, leads to the real cut function 

u B = G R (s,A(s, CO, C,C): 



•4$ = M w " ,,0 = «- ; YCkwiS(c,c) + v^gwi^(c,c) 

1 9zt 19 1 



(209) 



Later the linear versions are extensively used: 



•42 = y24 fl)scit) = . s - 



%/2 



4Wi(C,C) + V2^( s )y4.(c,c) 



r = «+i(^rzWi?i(c,c)-v^e?'(«)^(c,oV 



(210) 
(211) 
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5 Simple Applications 

In this section we give four simple examples of the use of shear-free and asymptotically shear- 
free NGCs in GR. The first is for asymptotically-linearized perturbations off the Schwarzschild 
metric, while the next two are from the class of algebraically-special metrics, namely the Robinson- 
Trautman metric and the type II twisting metrics; the fourth is for asymptotically static/stationary 
metrics. 



5.1 Linearized off Schwarzschild 

As a first example, we describe how the shear-free NGCs are applied in linear perturbations off the 
Schwarzschild metric. The ideas used here are intended to clarify the more complicated issues in the 
full nonlinear asymptotic theory. We will see that these linear perturbations greatly resemble our 
results from the previous section on the determination of the intrinsic center of charge in Maxwell 
theory, when there were small deviations from the Coulomb field. 

We begin with the Schwarzschild spacetime, treating the Schwarzschild mass, Ms c h = Mb, as 
a zeroth-order quantity, and integrate the linearized Bianchi identities for the linear Weyl tensor 
corrections. Though we could go on and find the linearized connection and metric, we stop just 
with the Weyl tensor. The radial behavior is given by the peeling theorem, so that we can start 
with the linearized asymptotic Bianchi identities, Equations (|56|) - ([58l) . 

Our main variables for the investigation are the asymptotic Weyl tensor components and the 
Bondi shear, a , with their related differential equations, i.e., the asymptotic Bianchi identities, 
Eq. (|56[) . ([57]) and (f55j). Assuming the gravitational radiation is weak, we treat cr° and &° as small. 
Keeping only linear terms in the Bianchi identities, the equations for ip® and ^ (the mass aspect) 
become 

$ + 3* = SV\ (212) 

* = 0, (213) 
f = f, (214) 

* = ^ + d 2 a°. (215) 

The ipi is small (first order), while the 

ir = if + + * y y 2 y + • • • ( 216 ) 
has the zeroth-order Schwarzschild mass plus first-order terms 

*° = -^-M Sch + 5*°, (217) 
fif 1 

¥ = -P\ (218) 

c 6 

In linear theory, the complex (mass) dipole moment, 

^Cferav.) = ^(mass) + J ' ( 219 ) 

is given |59j . on a particular Bondi cut with a Bondi tetrad (up to dimensional constants), by the 
I = 1 harmonic components of ip®, i.e., from the ip® 1 " in the expansion 

^=^%\ + ^i %% + ■■■ (220) 
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For a different cut and different tetrad, one needs the transformation law to the new and new 
i/j* oi . Under the tetrad transformation (a null rotation around n a ) to the asymptotically shear-free 
vector field, l* a , Equation (Tig) , 

r -> r a = l a - -m a - -m a + o(r- 2 ) } 

r r 

with, from Equations (|192j) and (|193|) . 

u B = C(r)L(CX)+C J (r)Y 2 %(CX) + --- (221) 
= ^C°(r) - ifW^CCC) + f J (r)^.(C, C) + • ■ • (222) 

i(«B,C,C) = fM^(C,C)-6f J (r)r4.(c,0 + ... (223) 

the linearized transformation is given by [6] 

= - (224) 

The extraction of the I = 1 part of ip®* should, in principle, be taken on the new cut given by the 
real u-q obtained from uq — £ a (r)? a (£, () + £ u (t)Y^ ; (C, C) + • • ■ with constant s in the expression, 
Equation (19"5|) . r = s + zA(s,£, £). However, because of the linearization, the extraction can be 
taken on the ub constant cuts. Following the same line of reasoning that led to the definition of 
center of charge, we demand the vanishing of the I = 1 part of . 
This leads immediately to 

$|i=i = 3L*| J= i, (225) 
or, using the decomposition into real and imaginary parts, ip1 z = ipin+iipi} and £ l (it rc t) = £i?( u rct) + 



4>IR = -^-M Sch e R (Uret), (226) 
= -^3^Msch£(Uret)- (227) 

Identifying |59i 137] the (intrinsic) angular momentum, either from the conventional linear iden- 
tification or from the Kerr metric, as 

r = S* = M Sch c£} (228) 

and the mass dipole as 

D l (mass} = Mb^Ch, (229) 
we have _ _ 

4? = -^^C( gr av) = -^{D\ mass) + ic _1 J*). (230) 

By inserting Equation (|230l) into Equation (|212[) . taking, respectively, the real and imaginary parts, 
using Equation (|2 18[) and the reality of "J, we find 

P l = M Sch CA = M Sch v l R , (231) 
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the kinematic expression of linear momentum and 

J" = 0, (232) 

the conservation of angular momentum. 

Finally, from the I = (0,1) parts of Equation (j2 1 3[) . we have, at this approximation, that the 
mass and linear momentum remain constant, i.e., M = M$ c b = Mb and S^ a = 0. Thus, we obtain 
the trivial equations of motion for the center of mass, 

Msch4" = 0- (233) 

It was the linearization that let to such simplifications. In Section |6j when nonlinear terms are 
included (in similar calculations), much more interesting and surprising physical results are found. 

5.2 The Robinson— Trautman metrics 

The algebraically-special type II Robinson-Trautman (RT) metrics are expressed in conventional 
RT coordinates, (r, r, f, Q, r now real, by [55] 

ds> =2(K- V Y r* + ^) dr 2 + 2drdr - (234) 

with 

K = 2V 2 P^d ( logVP , (235) 

Po = 1 + CC, (236) 

^2° = ^§(r). (237) 

The unknowns are the Weyl component ip^ (closely related to the Bondi mass), which is a function 
only of (real) r and the variable, V(t, £, C)j both of which satisfy the RT equation. (See below.) 
There remains the freedom 

t^t*= 9 (t), (238) 

which often is chosen so that V^ 7 ") = constant. However, we make a different choice. In the 
spherical harmonic expansion of V, 

V = v a l(C,() + v l iY 2 % + ..., (239) 

the r is chosen by normalizing the four- vector, v a , to one, i.e., v a v a — 1. The final field equation, 
the RT equation, is 

A' ~ 3^ ~ V 3 (S? T) BJ T) V - V- l t (T) V ■ df T) v) = 0. (240) 

These spacetimes, via the Goldberg-Sachs theorem, possess a degenerate shear-free PNV field, 
l a 7 that is surface- forming, (i.e., twist free). Using the tetrad constructed from l a we have that the 
Weyl components are of the form 

-00 = ipi = 0, 

4>2 ± o. 
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Furthermore, the metric contains a "real timelike world line, x a — £ a (i")," with normalized velocity 
vector v a — . All of these properties allow us to identify the RT metrics as being analogous to 
the real Lienard-Wiechert solutions of the Maxwell equations. 

Assuming for the moment that we have integrated the RT equation and know V = V(t, £, £), 
then, by the integral 

u = Jv(T,(X)dT = X RT (TX,0, (241) 

the UCF for the RT metrics has been found. The freedom of adding a(£, £) to the integral is just the 
supertranslation freedom in the choice of a Bondi coordinate system. From Xb,t(t, C>C) a variety 
of information can be obtained: the Bondi shear, er°, is given parametrically by 

a°( UB ,C,C) = gf T) X RT (r,C,C), (242) 

MB = Xr T (t,C,C), 

as well as the angle field L by 



L(u B ,(X) = 3 (t) Xrt(t,C,C), (243) 

In turn, from this information the RT metric (in the neighborhood of J + ) can, in principle, be re- 
expressed in terms of the Bondi coordinate system. In practice one must revert to approximations. 
These approximate calculations lead, via the Bondi mass aspect evolution equation, to both Bondi 
mass loss and to equations of motion for the world line, x a = ^ a (r). An alternate approximation for 
the mass loss and equations of motion is to insert the spherical harmonic expansion of V into the 
RT equation and look at the lowest harmonic terms. We omit further details aside from mentioning 
that we come back to these calculations in a more general context in Section [6] 

5.3 Type II twisting metrics 

It was pointed out in the previous section that the RT metrics are the general relativistic analogues 
of the (real) Lienard-Wiechert Maxwell fields. The type II algebraically-special twisting metrics 
are the gravitational analogues of the complex Lienard-Wiechert Maxwell fields described earlier. 
Unfortunately they are far more complicated than the RT metrics. In spite of the large literature 
and much effort there are very few known solutions and much still to be learned [551 EH 122] • We 
give a very brief description of them, emphasizing only the items of relevance to us. 

A null tetrad system (and null geodesic coordinates) can be adopted for the type II metrics so 
that the Weyl tetrad components are such that 

■00 = i/ji = 0, 

02 + 0. 

It follows from the Goldberg-Sachs theorem that the degenerate principal null congruence is geodesic 
and shear-free. Thus, from the earlier discussions it follows that there is a unique angle field, 
L(ub,CiO- As with the complex Lienard-Wiechert Maxwell fields, the type II metrics and Weyl 
tensors are given in terms of the angle field, L(mb,CiO- I n fact, the entire metric and the field 
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equations (the asymptotic Bianchi identities) can be written in terms of L and a Weyl tensor 
component (essentially the Bondi mass). Since L(uqXX) describes a unique shear-free NGC, it 
can be written, parametrically, in terms of a unique GCF, namely a UCF, -X"( type ii)( t j Cj C): i- e -> we 
have that 

L(UB,(X) = 9(T)^(typc II), 

«B = -^(typc II) ( r ; C, 0- 

Since X( typc m(r, £, C) can be expanded in spherical harmonics, the I = (0,1) harmonics can be 
identified with a (unique) complex world line in "H-space. The asymptotic Bianchi identities then 
yield both kinematic equations (for angular momentum and the Bondi linear momentum) and equa- 
tions of motion for the world line, analogous to those obtained for the Schwarzschild perturbation 
and the RT metrics. As a kinematic example, the imaginary part of the world line is identified as 
the intrinsic spin, the same identification as in the Kerr metric, 

S 4 = MbC#. (244) 

In Section [6l a version of these results will be derived in a far more general context. 



5.4 Asymptotically static and stationary spacetimes 

By defining asymptotically static or stationary spacetimes as those asymptotically-flat spacetimes 
where the asymptotic variables are 'time' independent, i.e., «b independent, we can look at our 
procedure for transforming to the complex center of mass (or complex center of charge). This 
example, though very special, has the huge advantage in that it can be done exactly, without the 
use of perturbations [2] • 

Imposing time independence on the asymptotic Bianchi identities, Equations (|56p - (|58| . 



4° = -Bl^ + Aft (245) 

i'i = -3^2+2^3°, (246) 

= _g^o + 3{r o^ j (247) 

and reality condition 

* = V>2 + + crW = ^, (248) 
we have, using Equations (fBTj) and (|52|) with a — 0, that 



^3° 


= V>4° = o, 


(249) 


3^2° 


= o, 


(250) 




= 3a>°, 


(251) 




= ^+3 2 a = V^ + 3 2 ( r = '*. 


(252) 



From Equation (|252[) . we find (after a simple calculation) that the imaginary part of i/j® * s 
determined by the 'magnetic' [44] part of the Bondi shear (spin-weight s = 2) and thus must 
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contain harmonics only of I > 2. But from Equation (I250p . we find that ip® contains only the / = 
harmonic. From this it follows that the 'magnetic' part of the shear must vanish. The remaining 
part of the shear, i.e., the 'electric' part, which by assumption is time independent, can be made to 
vanish by a supertranslation, via the Sachs theorem: 

u B = u B +a((,(), (253) 

5(CC) = a(C,C)+3 2 a(C,C). 

In this Bondi frame, (i.e., frame with a vanishing shear), Equation (|25ip . implies that 

= _ (254) 

V'l 1 = ^2 -^C(grav) = ^2 (^(mass) + ic (255) 

using the conventionally accepted physical identification of the complex gravitational dipole. (Since 
the shear vanishes, this agrees with probably all the various attempted identifications.) 
From the mass identification, becomes 

o = _^ Mb . (256) 

Since the Bondi shear is zero, the asymptotically shear-free congruences are determined by the same 
GCFs as in flat spaces, i.e., we have 

L(u B ,(X) = 3 (r) G(T,C,C)=C(T)m Q (C,a (257) 

u B = C(r)i a (CX)- (258) 

Our procedure for the identification of the complex center of mass, namely setting ipl° = in 
the transformation, Equation (|276[) , 

leads, after using Equations (|254|) . (|249|) and f)257f) . to 

V>? = 3LV 2 °, (259) 
_ 6V2G 

- c 2 ^C(grav)' 

^C(grav) = MbC- 

From the time independence, £ l , the spatial part of the world line is a constant vector. By a 
(real) spatial Poincare transformation (from the BMS group) , the real part of £ l can be made to 
vanish, while by ordinary rotation the imaginary part of can be made to point in the three- 
direction. Using the the gauge freedom in the choice of r we set £°(r) = r. Then pulling all 
these items together, we have for the complex world line, the UCF, L(ub,(X) an d the angular 
momentum, J 1 : 

C(t) = (r, 0,0,1^), (260) 
u B = C(r)h 



UB = C(T)fa(C,C)^-^ 3 n°3, 



J 1 = S t = M B ^Sl=M B c(0,0,t) = M B (ti- 
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Thus, we have the complex center of mass on the complex world line, z a = £ a (r). 

These results for the lower multipole moments, i.e., I = 0,1, are identical to those of the Kerr 
metric. The higher moments are still present (appearing in higher r _1 terms in the Weyl tensor) 
and are not affected by these results. 



6 MAIN RESULTS 



51 



6 Main Results 

We saw in Sections [3] and [4] how shear-free and asymptotically shear-free NGCs determine arbitrary 
complex analytic world lines in the auxiliary complex 'H-space (or complex Minkowski space). In 
the examples from Sections |3] and El we saw how, in each of the cases, one could pick out a 
special GCF, referred to as the UCF, and the associated complex world line by a transformation 
to the complex center of mass or charge by requiring that the complex dipoles vanish. In the 
present section we consider the same problem, but now perturbatively for the general situation 
of asymptotically-flat spacetimes satisfying either the vacuum Einstein or the Einstein-Maxwell 
equations in the neighborhood of future null infinity. Since the calculations are relatively long and 
complicated, we give the basic details only for the vacuum case, but then present the final results 
for the Einstein-Maxwell case without an argument. 

We begin with the Reissner-Nordstrdm metric, considering both the mass and the charge as 
zeroth-order quantities, and perturb from it. The perturbation data is considered to be first order 
and the perturbations themselves are general in the class of analytic asymptotically-flat spacetimes. 
Though our considerations are for arbitrary mass and charge distributions in the interior, we look at 
the fields in the neighborhood of 3 + . The calculations are carried to second order in the perturbation 
data. Throughout we use expansions in spherical harmonics and their tensor harmonic versions, 
but terminate the expansions after 1 = 2. Clebsch-Gordon expansions are frequently used. See 
Appendix [C] 



6.1 A brief summary — Before continuing 

Very briefly, for the purpose of organizing the many strands so far developed, we summarize our 
procedure for finding the complex center of mass. We begin with the gravitational radiation data, 
the Bondi shear, er°(uB, £, C) an d solve the good-cut equation, 

<5 2 Z = a°(ZXX), 

with solution ub = Z(z a and the four complex parameters z a defining the solution space. 

Next we consider an arbitrary complex world line in the solution space, z a = £ a (r) = (£°(t), 

so that Mb = Z{^ a (r) : £, = G(r, £, (), a GCF, which can be expanded in spherical harmonics as 

u B = G(TX,C)=C(r)ia<:X + e'(r)Y 2 %+... (261) 



y/2 2 

Assuming slow motion and the gauge condition £°(t) = r (see Section QJ, we have 

r 1 
71 ~ 2 



u B = — -C (r)Y° + C 3 (r)Y 2 % + ... (262) 



(Though the world line is arbitrary, the quadrupole term, £ y (7~), and higher harmonics, are deter- 
mined by both the Bondi shear and the world line.) 
The inverse function, 

t = T(u xet ,C,0, 

M rct = V2u Bl 
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can be found by the following iteration process [37]; writing Equation (|262[) as 

T = u Tet + F(r,C,C), (263) 

with 

F(r, C, C) = ^ (t)YS (C, - V2C j (r)Y? l3 (C, C) + • • • , (264) 
the iteration relationship, with the zeroth-order iterate, To = u rc t, is 

+ F(r„_ 1 ,C,C)- (265) 

Though the second iterate easily becomes 

r = T(u IcU C, C) = u Iet + F (u rot + F{u rct XX), C, C) « "rot + i 7 + ^d„ r6t F. (266) 
For most of our calculations, all that is needed is the first iterate, given by 

2 



T(u iet , C, C) = "ret + ^rC{u iet )Y°(C, - v / 2f J '( Wret )^ i (C, 0- (267) 



This relationship is important later. 

We also have the linearized reality relations - easily found earlier or from Equation (I267p : 

t = s + i\, (268) 
A = A(sX,0 = ^}(s)Y 1 i -V2^(s)Y 2 %, (269) 

r = s + i(^$(s)Y°-V2g(s)Yl\ (270) 

2 

The associated angle field, L, and the Bondi shear, a , are given parametrically by 



^scR) = V2G fl ( S ,C,C) = V^««^^^ (2n) 



L( UB ,C,C) = S (T )G(r,C,C) (272) 
= e(r)Y^-6C j (r)Y^ + ... 

and 



c7 (r,C,C) = 3^ T )G(r,C,C), (273) 
= 24^(^4-+..., 

using mb = G, while the asymptotically shear-free NGC is given (again) by the null rotation 



l* a = l a + bm a + bm a + bbn a , (274) 

m* a = m a + bn a , 

n* a = n a , 

b = -L/r + 0(r- 2 ). 
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r ° 


= 


4LVi*° 


+ 6L 2 ?p* a 




= 


3L<ip% -+ 


- 3i 2 V-3 - 




= v> 2 °- 


2L^3 4 






= ^l- 






v> 4 ° 


= 







From Equation (|274[) . the transformed asymptotic Weyl tensor becomes, Equations (|275[) - (|279[) . 

4L^f + L 4 ^, (275) 

(276) 
(277) 
(278) 
(279) 

The procedure is centered on Equation (|276j) . where we search for and set to zero the I = 1 
harmonic in on an s = constant slice. This determines the complex center-of-mass world line 
and singles out a particular GCF referred to as the UCF, 

X(t,(,0 = G(t,(,0, 

with the real version, 

X R (s,£,Q = G R (s,£,0, (280) 

for the gravitational field in the general asymptotically-flat case. 

(For the case of the Einstein-Maxwell fields, in general, there will be two complex world lines, 
one for the center of charge, the other for the center of mass and the two associated UCFs. For 
later use we note that the gravitational world line will be denoted by £ a , while the electromagnetic 
world line by r\ a . Later we consider the special case when the two world lines and the two UCFs 
coincide, i.e., £ a = if.) 

From the assumption that er° and L are first order and, from Equation (|52[) . that ip® = 3(7 , 
Equation (|276j) . to second order, is 

= ^-3L[*-5 2 ct°], (281) 

where ip^ has been replaced by the mass aspect, Equation (|54|) . "J w tp® + 3 2 it . 
Using the spherical harmonic expansions (see Equations (|272j) and (|273p ). 

* = *° + + <S y Y 2 °« + • • • , (282) 

& = ^ 1 %\+^ ij Yi ij + ..., (283) 

r * = ri°%\ + ^ 0ii Yi ij + ..., (284) 

L(u B X,0 = e(r)Y^ i -6e j (r)Y 2 1 ij +..., (285) 

a°(u B ,C,0 = 24f'(r)y a ^+... (286) 

Remembering that "J -0 is zeroth order, Equation ((281]), becomes 



2»j 

-3[f(r)>H 6e\r)Y 2 \^° + **** + - 24|* J (r)}l^.j 



or 
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#° = r^ + ^Y^-ZeirW^-W^iT^YV (287) 
-3e(rW> -24f(r)]Y 1 \Y 2 % 

+18C 3 (r)Y 2 \^° + * k lKHr)Y 2 \ 3 Y? k + 18£ w (r)[*« - 24? j (T)}Yi kl Y 2 % . 
Note that the right-hand side of Equation (|287|) depends initially on both r and u m %, with r = 

This equation, though complicated and unattractive, is our main source of information concern- 
ing the complex center-of-mass world line. Extracting this information, i.e., determining -0^°* at 
constant values of s by expressing r and it rc t as functions of s and setting it equal to zero, takes 
considerable effort. 

6.2 The complex center-of-mass world line 

Before trying to determine V^ 01 , several comments are in order: 

1. As mentioned earlier, the right-hand side of Equation (|287p is a function of both t, via the 
the GC variables, £*, £ u and the u rc t, via the 4>?, ^x 1 an d V?. The extraction of the I = 1 
part of ip*° must be taken on the constant V cuts. In other words, both r and u rc t must be 
eliminated by using Equations P58" ]) - (|2TT ]) and ([557]) . 

2. This elimination must be done in the linear terms, e.g., from Equation (|271l) . 

« »?(«) - - 2^( s )r 2 ^ 

or, from Equation (|270[) , 
In the nonlinear terms we can simply use 

""rot = T = S. 

3. In the Clebsch-Gordon expansions of the harmonic products, though we need both the I = 1 
and I = 2 terms in the calculation, we keep at the end only the I = 1 terms for the ip® 1 . (Note 
that there are no I — terms since ipi° is spin weight s = 1.) 

4. The calculation to determine the for the constant V slices was probably the most tedious 
and lengthy in this work. The importance of the results necessitated the calculation be 
repeated several times. 
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Expanding and organizing Equation (|287p with the linear terms given explicitly and the quadratic 
terms collected in the expression for A, we obtain the long expression with all terms functions of 
either r or u rot : 



+ v ( l ij Y> - 3*°^ + 18*°^, + A 



A 



(288) 
(289) 



kl 



24^(r) 



yl yO 
I 2ij I 2kl 



= A'Y^ + A^. 

Using Equations (|269[) and (|271[) . the u TCt and r in Equation (|288|) are replaced by s. On the 
right-hand side all the variables, e.g., if)® , £ fe , etc., are functions of 's'; their functional forms are 
the same as when they were functions of u rct and r; the linear terms are again explicitly given and 
the quadratic terms are collected in the Bjj + Bt + A 



CO 

Bu 



+V2^'$Y 2 ° kl Y 2 \ 



(290) 



2ij: 



B T = z3 V2*o ( _ W (s)y o yi. + (^yi + 



2 (J 



ay 

To proceed, we use the complex center-of-mass condition, namely, ipl l>k (s) — 0, and solve for 
i/ji(u Ie t)- This is accomplished by first reversing the calculation via 

and then, before extracting the I = 1 harmonic component, replacing the s by u re t, via the inverse 
of Equation ([2TTj) . 



,((iJ)scfl) 



o SiiV"ret 



2 -s«\-ret Hli V^^^Urct ) I 2ij> 



using 



(291) 



X R (s,C,C)=Gii( S ,C,C)- 
In this process several of the quadratic terms cancel out and new ones arise. 

The final expression for ip^ , given in terms of the complex world line f expressed as a function 
of u re t, then becomes our basic equation: 

3y/2 



Oj _ 



■fig 



18 



= 3*T-y£V 2 



108 



■5 



(292) 



3^2 ^ 



24 



216V2 



6 MAIN RESULTS 



56 



This, which becomes the analogue of the Newtonian dipole expression D = MR, is our central 
relationship. Almost all of our results in the following sections follow directly from it. 

We emphasize that prior to this discussion/derivation, the ip^ and the £ J were independent 
quantities but in the final expression the ip^ is now a function of the f' . 

Note that the linear term 

ip? = 3*°f 

coincides with the earlier results in the stationary case, Equation (|259|) . From 



we have 



^ = _^ MbCj (293) 



6V2G 



Dl 



c 2 ^C(grav) 



-^(grav) = M b(Cr +i&)= (£>(mass) + ic J*)- 

We will see shortly that there is a great deal of physical content to be found in the nonlinear terms 
in Equation ([292]) . 



6.3 Results 
6.3.1 Preliminaries 

Before describing in details our results, it appears to be worthwhile to very roughly survey the 
results and describe the logical steps taken to reach them. Virtually everything, as we said earlier, 
follows from the equation for tp^ 1 , i.e., Equation (I292p . 

In the final results, we will include, from Section |3l the Maxwell field. Though we do consider 
the case where the two complex world lines (the complex center-of-mass and center-of-charge lines) 
differ from each other, some discussion will be directed to the special case of coinciding world lines. 

1. The first step is to decompose the ip® 1 into its real and imaginary parts, identifying the real 
center of mass and the total angular momentum (as seen from infinity) described in the 
given Bondi coordinate system. In a different Bondi system they would undergo a specific 
transformation. These results are the analogues of Equation (|II5[) . It should be emphasized 
that there are alternative definitions, |59j . but using our approximations all should reduce to 
our expression. 

2. The second step is to look at the evolution equation for ip® 1 , i.e., insert our ip® into the 
Bianchi identity, (|57|). now including the Maxwell field. We obtain, from the I — 1 harmonic, 
the evolution of ip^ 1 . After again decomposing it into the real and imaginary parts we find 
the kinematic description of the Bondi linear momentum, P l = M£ 1 ' + . . . (i.e., the usual 
kinematic expression P = Mv plus additional terms) and the evolution (conservation law) for 
the angular momentum including a flux expression, i.e., J 1 ' = (Flux) 1 . 
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3. The third step is to reinsert the kinematic expression for the Bondi mass into the evolution 
equation for ip®, i.e., Equation (1561) . From the reality condition on the mass aspect, Equa- 
tion (|55l) or Equation (|53l) . only the real part is relevant. It leads to the evolution equation 
for the real part of the complex world line, a second-order ODE, that can be identified with 
Newton's second law, F l — M^ 1 " , with F l being the recoil and radiation reaction forces. The 
I = harmonic term is the energy /mass loss equation of Bondi. 

Before continuing we note that the I = 2 coefficients in i/j® and tp®, i.e., and ip 2 l * \ appear 
frequently in second-order expressions, e.g., in Equation (|292[) . Thus, knowing them, in terms of 
the free data, to first order is sufficient. By going to the linearized Bianchi identities (with the 
linearized Maxwell field) and the expression for the Bondi shear a , 

$ = -oV, (294) 

^ = -&l% + 2kq%, (295) 

<7°(r,C,0 = 24^(t)Y 2 %+..., (296) 

$ = g, (297) 



1° _ r>„^ill/„, \ V 1 ±p i i'"vl 

3 



, a - -2qrr"{u Tet )Y^--Qi Y 2 \ p (298) 



we easily find 



W = -24C(u Iet ), (299) 

where a constant of integration was set to zero via initial conditions. These expressions are fre- 
quently used in the following. In future expressions we will restore explicitly 'c', via the derivative, 
(') — > c _1 (') and replace the gravitational coupling constant by k — 2Gc -4 . 

6.3.2 The real center of mass and the angular momentum 

Returning to the basic relation, Equation (|292|) . using Equations (|299p we obtain 

i& = 3*°e + &Lf j e + ^c-^vv (300) 

5 5 
— g-* 4 h g — + — g— c (j<iQc 5i? 

+i ^3a^V2 r , r _ | Gc -7^'"^ £ifci . 

By replacing "J and 'i' 2 , in terms of the Bondi mass and linear momentum, then decomposing the 
individual terms, e.g., £ l = + = Q 1 ^ + iQ 1 ^, into their real and imaginary parts, the full 
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expression is decomposed as 

(mass) 



^ = -^-(DL^+ic-iT) (301) 



= - ( ^(M B & + iM B %) + ... 
cr 

The physical identifications - first from the real part, are, initially, a tentative definition of the mass 
dipole moment, 

D { Si) = ^B^ + c-^MB^ + iMB^j^+P^, (302) 



5 



36V2 ! 2 j jj \/2 5 i ij„, 

J 2(12) ^—12 tim ^rnfc ^ „— 5 „cmk /-\im ill \ „ 

- s — g — G c ?/ Si? - q£n Qm > € ik] 

36-\/2 _ 4 2 , -4„2>i fe// 

and - from the imaginary part, again, a tentative definition of the total angular momentum, 

J( T ^" = M B c# + {m b &'& + ^M B &f} + (304) 

J* j - - 3 -^^^mt-^)-^M B ^f (305) 
6V2 , , t jiMr , „-&„£% n ij»> & „-i„£mkriimin 



c- 3 q 2 V t R , '^+c- 3 q 2 e R V R / 'e lkr 



The reason for referring to these identifications as tentative is the following: 
If there were no Maxwell field present, then the terms involving the electromagnetic dipole, qrf , 
and quadruple, Q™, would not appear and these identifications, D^j^s and J^ 1 , would then be 
considered to be firm; however, if a Maxwell field is present, we will see later that the identifications 
must be modified. Extra Maxwell terms are 'automatically' added to the above expressions when 
the conservation laws are considered. 

As an important point we must mention a short cut that we have already taken in the interests of 
simplifying the presentation. When the linearized Equation (|300[) is substituted into the linearized 
Bianchi identity, ip® = — cS^S?, we obtain the linear expression for the momentum, P l = — ^^2% m 
terms of the linearized expression for ipi 1 , namely, 

P l = M B $ - \c-\*nt = Mb^ - f <T 3 «. (306) 
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This expression is then 'fed' into the full nonlinear Equation (|300p leading to the relations Equa- 
tions (I3U2T) and ([304]) . 

Considering now only the pure gravitational case, there are several comments and observations 
to be made. 

1. Equations fl302[) and (I304|) have been split into two types of terms: terms that contain only 
dipole information and terms that contain quadrupole information. The dipole terms are 
explicitly given, while the quadrupole terms are hidden in the T>& 1 and 1 . 

2. In J 1 we identify S l as the intrinsic or spin angular momentum, 

S i = M B c£l . (307) 
This identification comes from the Kerr or Kerr-Newman metric [37] ■ The second term, 

MB&'&eay = x ify, (308) 
is the orbital angular momentum. The third term, 

though very small, represents a spin-spin contribution to the total angular momentum. 

(T) i 

3. In the mass dipole expression D( mass y the first term is the classical Newton mass dipole, while 
the next two are dynamical spin contributions. 

In the following Sections 16.3.31 and 16.3.41 further physical results (with more comments and 
observations) will be found from the dynamic equations (asymptotic Bianchi identities) when they 
are applied to and j( T )\ 

6.3.3 The evolution of the complex center of mass 

The evolution of the mass dipole and the angular momentum, defined from the tp^, Equation (|301[) . 
is determined by the Bianchi identity, 

i>l = -902 + 2frVa + 2fc0°0g. (309) 

In the analysis of this relationship, the asymptotic Maxwell equations 

4>l = -8$, (310) 

and their solution, from Section |3j Equation (|123[) (needed only to first order), 

$ = 2qr ) i (u Tet )Y 1 1 i +c- 1 Q i i%\ j + ..., (311) 
0? = q + V2c- 1 qr ] i '(u Iet )Y 1 i + ^-c- 2 Qf'Y° j + ..., 
$ = -2c-W>re t )^7 1 -^- 3 Q«"Y 2 Ti + ..., 
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must be used. By extracting the I — 1 harmonic from Equation (|309p . we find 

^ ' = V2cV + ?2(12 f^ ^f J 'e w - SVlGc-Vv 1 " + l iV2Gc- e q 2 ^ "r, m >e mji (312) 





Using Equation pOlj) . 
with the (real) 



MB -^" (313) 

P* = (314) 

we obtain, (1) from the real part, the kinematic expression for the (real) linear momentum and, (2) 
from the imaginary part, the conservation or flux law for angular momentum. 
(1) Linear Momentum: 

P l = D^-^c-tq^' + W, (315) 

IT = - ^fH, - ^{«' - 

2v^2 — 5„(„kn / ~ l kiii , m kii r^kin im r^im ill , -ml nimlll^ 



15 



+ «X - vrQT + viQm"} 



Using 



I ■* n m i " n ml1 " n m i"n ml '"\ c 
45 "^a/ ^ B ^ £ ^ M ' l i l 



DfJ: s ' s) = - c- 1 {m b #£ + i.W,^;'}' e jki + V 



(mass) > 



we get the kinematic expression for the linear momentum, 



P l = M B Zk'-lc-Vv%-c-HMBZ%&+M B &"8 (316) 

+^ C - 5 <z{^"Qr + r?fQS" - Vr'Qe 1 ''" + vT'Qm'"} 

I "* _— e r n m ^ " n ml1 " n m i " n mllll \ c 
45 ^ B ^Af / e 0'i' 
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J* = M B £A' " |c-V»# + c- 1 M B ^ , c iifc +S i ) (317) 

s i = -c-H^B^+^M B (^'y}ei«+^5: a ' ss) +ir, 

(2) Angular Momentum Flux: 

J l ' = {Flux)\ (318) 
J' = J^ + ^c" 2 ^', (319) 

J^ 1 = M B cei + {M B ^U+^Bekf}^i + J #i , (320) 

(Fluxf = -lM!g{^^''+^^ 7 }e yi (321) 

2\/2 _ At hii^uin kiir^kiii 



15 
2a/2 



<jc *{r,?'QT'" ~ VrQm'"}- 



15 

There are a variety of comments to be made about the physical content contained in these 
ations: 



• The first term of P l is the standard Newtonian kinematic expression for the linear momentum, 

• The second term, — ^c~ 3 q 2 ri R ' , which is a contribution from the second derivative of the 
electric dipole moment, qrf R , plays a special role for the case when the complex center of mass 
coincides with the complex center of charge, rf = £ a . In this case, the second term is exactly 
the contribution to the momentum that yields the classical radiation reaction force of classical 
electrodynamics [30] . 

• The third term, c" 1 Mb^^j 'tijk, is the classical Mathisson-Papapetrou spin- velocity contri- 
bution to the linear momentum. If the evolution equation (angular momentum conservation) 
is used, this term becomes third order. 

• Many of the remaining terms in P 1 , though apparently second order, are really of higher order 
when the dynamics are considered. Others involve quadrupole interactions, which contain high 
powers of c _1 . Though it is nice to see the Mathisson-Papapetrou term, it should be treated 
more as a suggestive result rather than a physical prediction. 

Aside: Later, in the discussion of the Bondi energy loss theorem, we will see that we can 
relate £ lJ , i.e., the I = 2 shear term, to the gravitational quadrupole by 

e j = «g + <#) = ^f^iQZL + <L) = ^ Gc 2 ^ ;iv - (322) 
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• In the expression for J 1 we have already identified, in the earlier discussion, the first two terms, 
Mbc£j and Mb^' CR e ikj as the intrinsic spin angular momentum and the orbital angular 
momentum. The third term, ^M^}^' e^-j — 2 m b c 1 e ikji a spin-spin interaction term, 
considerably smaller, can be interpreted as a spin-precession contribution to the total angular 
momentum. An interesting contribution to the total angular momentum comes from the term, 
^c~ 2 q 2 n}' = ^c~ 2 qD\' { , i.e., a contribution from the time-varying magnetic dipole. 

• As mentioned earlier, our identification of j( T ' 1 as the total angular momentum in the ab- 
sence of a Maxwell field agrees with most other identifications (assuming our approxima- 
tions). Very strong support of this view, with the Maxwell terms added in, comes from the 
flux law. In Equation Q321I) we see that there are five flux terms, the first is from the grav- 
itational quadrupole flux, the second and third are from the classical electromagnetic dipole 
and electromagnetic quadrupole flux, while the fourth and fifth come from electromagnetic- 
gravitational coupling. The Maxwell dipole part is identical to that derived from pure Maxwell 
theory [30j . We emphasize that this angular momentum flux law has little to do directly 
with the chosen definition of angular momentum. The imaginary part of the Bianchi identity, 
Equation (|309l) . is the conservation law. How to identify the different terms, i.e., identifying 
the time derivative of the angular momentum and the flux terms, comes from different argu- 
ments. The identification of the Maxwell contribution to total angular momentum and the 
flux contain certain arbitrary assignments: some terms on the left-hand side of the equation, 
i.e., terms with a time derivative, could have been moved onto the right-hand side and been 
called 'flux' terms. However, our assignments were governed by the question of what terms 
appeared most naturally to be explicit time derivatives, thereby being assigned to the time 
derivative of the angular momentum. 

• The angular momentum conservation law can be considered as the evolution equation for the 
imaginary part of the complex world line, i.e., £}(u rc t)- The evolution for the real part is 
found from the Bondi energy-momentum loss equation. 

• In the special case where the complex centers of mass and charge coincide, n a — £ a , we have a 
rather attractive identification: since now the magnetic dipole moment is given by D l M = q£} 
and the spin by S l = M^cQ, we have that the gyromagnetic ratio is 

\S l \ = Mbc 

leading to the Dirac value of <?, i.e., g = 2. 

6.3.4 The evolution of the Bondi energy- momentum 

Finally, to obtain the equations of motion, we substitute the kinematic 
Bondi evolution equation, the Bianchi identity, Equation (|56[) : 

4° = -g^ 3 ° + CT V2 + 

or its much more useful and attractive (real) equivalent expression 



expression for P % into the 
(323) 

(324) 
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with 

# = V" + S 2 ^ + <fU = *. (325) 
Remark: The Bondi mass, Mb = ~~ 2 \/2G ' anc ^ ^ e or i§ ma l mass of the Reissner-Nordstrdm 
(Schwarzschild) unperturbed metric, Mrn = — 9 ^ G ^2 ° 1 i- e -' the / = harmonic of ipQ, differ by 

a quadratic term in the shear, the / = part of era. This suggests that the observed mass of an 
object is partially determined by its time-dependent quadrupole moment. 

Looking only at the I — and I — 1 spherical harmonics and switching to the "' derivatives with 
the c _1 inserted, we first obtain the Bondi mass loss theorem: 



Mi 2(12) c / ij ij jj j,\ 2q 1 in in , in iu\ f ( r^ii ' 1 " 1 rV>i "> ' r\i] "' 

b = + «j ■ «j J - 3^ (»?5* Vr +Vi Vi )- j^g? [Qe Qe + Qm Qm 



5G ■ -1 j 3c 5 vi« '« ■ ■« / isoc 

If we identify £ lJ with the gravitational quadrupole moment QQ rav via 



1 j 11 



4 



and the electric and magnetic dipole moments by 

the mass loss theorem becomes 

Mb = ~5^7F (^Mass^Mass + QspLQspin) (326) 



5c 7 



M 



The mass/energy loss equation contains the classical energy loss due to electric and magnetic 
dipole radiation and electric and magnetic quadrupole {Q%,Q l ^j) radiation. The gravitational en- 
ergy loss is the conventional quadrupole loss by the above identification of £ y with the gravitational 
quadrupole moment QQ rav - 

The momentum loss equation, from the I = 1 part of Equation (|324[) . becomes 

P k ' = F lecoil , (327) 

F k _ 2G ( n lj "'n ij "' n lj/ " n ij "'\ r ■„ —n l 'r)i' P .„ 

r recoil — ^5 C 6 V S pin Mass ^Mass^Spin^/ tllk }iI E dk 



5c 5 y E E U M U M j 54 q c 6 ^ E M ^ 



15c 5 

Finally, substituting the P 1 from Equation (|317[) . we have Newton's second law of motion: 

M B &" = F*, (328) 

with 

F l = -M^H - c-iM^tf'Yem + \c-* q 2 r?>" + - 5". (329) 

There are several things to observe and comment on concerning Equations (|328l) and (13291) : 
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• If the complex world line associated with the Maxwell center of charge coincides with the 
complex center of mass, i.e., if rf = £ l , the term 



becomes the classical electrodynamic radiation reaction force. 

• This result follows directly from the Einstein-Maxwell equations. There was no model building 
other than requiring that the two complex world lines coincide. Furthermore, there was no 
mass renormalization; the mass was simply the conventional Bondi mass as seen at infinity. 
The problem of the runaway solutions, though not solved here, is converted to the stability 
of the Einstein-Maxwell equations with the 'coinciding' condition on the two world lines. If 
the two world lines do not coincide, i.e., the Maxwell world line forms independent data, then 
there is no problem of unstable behavior. This suggests a resolution to the problem of the 
unstable solutions: one should treat the source as a structured object, not a point, and centers 
of mass and charge as independent quantities. 

• The F* ccoil is the recoil force from momentum radiation. 

• The S 1 ' = — F^ft can be interpreted as the gravitational radiation reaction. 

• The first term in F l , i.e., — Afg£]/, is identical to a term in the classical Lorentz-Dirac equa- 
tions of motion. Again it is nice to see it appearing, but with the use of the mass loss equation 
it is in reality third order. 

6.4 Other related results 

The ideas involved in the identification, at future null infinity, of interior physical quantities that 
were developed in the proceeding sections can also be applied to a variety of different perturbation 
schemes. Bramson, Adamo and Newman OH [3] have investigated how gravitational perturbations 
originating solely from a Maxwell radiation field can be carried through again using the asymptotic 
Bianchi identities to obtain, in a different context, the same identifications: a complex center-of- 
mass/charge world line, energy and momentum loss, as well as an angular momentum flux law that 
agrees exactly with the predictions of classical electromagnetic field theory. This scheme yields (up 
to the order of the perturbation) an approximation for the metric in the interior of the perturbed 
spacetime. 

We briefly describe this procedure. One initially chooses as a background an exact solution 
of the Einstein equations; three cases were studied, flat Minkowski spacetime, the Schwarzschild 
spacetime with a 'small' mass and the Schwarzschild spacetime with a finite, 'zero order', mass. For 
such backgrounds, the set of spin coefficients is known and fixed. On this background the Maxwell 
equations were integrated to obtain the desired electromagnetic field that acts as the gravitational 
perturbation. Bramson has done this for a pure electric dipole solution [9l [T] on the Minkowski 
background. Recent work has used an electric and magnetic dipole field with a Coulomb charge [3] . 
The resulting Maxwell field, in each case, is then inserted into the asymptotic Bianchi identities, 
which, in turn, determine the behavior of the perturbed asymptotic Weyl tensor, i.e., the Maxwell 
field induces nontrivial changes to the gravitational field. Treating the Maxwell field as first order, 
the calculations were done to second order, as was done earlier in this review. 




i ill 



(330) 
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Using the just obtained Weyl tensor terms, one can proceed to the integration of the spin- 
coefficient equations and the second-order metric tensor. For example, one finds that the dipolc 
Maxwell field induces a second-order Bondi shear, a . (This in principle would lead to a fourth-order 
gravitational energy loss, which in our approximation is ignored.) 

Returning to the point of view of this section, the perturbed Weyl tensor can now be used to 
obtain the same physical identifications described earlier, i.e., by employing a null rotation to set 
ipi = 0, equations of motion and asymptotic physical quantities, (e.g., center of mass and charge, 
kinematic expressions for momentum and angular momentum, etc.) for the interior of the system 
could be found. Although we will not repeat these calculations here, we present a few of the results. 
Though the calculations are similar to the earlier ones, they differ in two ways: there is no first-order 
freely given Bondi shear and the perturbation term orders are different. 

For instance, the perturbations induced by a Coulomb charge and general electromagnetic dipole 
Maxwell field in a Schwarzschild background lead to energy, momentum, and angular momentum 
flux relations [3]: 




(331) 



P 



——n k "n j "e, 



J' 




all of which agree exactly with predictions from classical field theory [30] . 

The familiarity of such results is an exhibit in favor of the physical identification methods 
described in this review, i.e., they are a confirmation of the consistency of the identification scheme. 
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7 Gauge (BMS) Invariance 

The issue of gauge invariance, the understanding of which is not obvious or easy, must now be 
addressed. The claim is that the work described here is in fact gauge (or BMS) invariant. 

First of all we have, or its real part, 3 + . On 3^, for each choice of space-time interior 
and solution of the Einstein-Maxwell equations, we have its UCF, either in its complex version, 
ub = X(t, C, £), or its real version, Eq. (|291[) . The geometric picture of the UCF is a one-parameter 
family of slicings (complex or real) of 3 J or 3 + . This is a geometric construct that has a different 
appearance or description in different Bondi coordinate systems. It is this difference that we must 
investigate. We concentrate on the complex version. 

Under the action of the supertranslation, Eq. flBTl) . we have: 



X(r, C, -► X(r, C, C) = X(t, C, C) + a(C, 0, (332) 

with a(C, C) an arbitrary complex smooth function on (complexified) S 2 . Its effect is to add on a 
constant to each spherical harmonic coefficient. The special case of translations, with 

a(C,0=d a i a (C,0, (333) 

simply adds to the I = (0, 1) harmonic components the complex constants d a , so, via Eq. (|26ip . we 
have the (complex) Poincare translations, 

£" -> £ a = £° + d a . (334) 
The action of the homogeneous Lorentz transformations, Eq. (l68p . 

ub = Kub, (335) 

1 + cc 

K = _ _ ^ =— (336) 

(a( + b)(a( + b) + (c( + d){c( + d) 

C - ad-bc=l. (337) 

cC + d 

e x = <±d (338) 
cC + d 

is considerably more complicated. It leads to 

X(r,C,C)=^(r,C,0- (339) 

Before discussing the relevant effects of the Lorentz transformations on our considerations we 
first digress and describe an important technical issue concerning representation of the homogeneous 
Lorentz group. 

The representation theory of the Lorentz group, developed and described by Gelfand, Graev and 
Vilenkin [IB] used homogeneous functions of two complex variables (homogeneous of degrees, rii — 1 
and ri2 — 1) as the representation space. Here we summarize these ideas via an equivalent method 
[2T1 [15] using spin- weighted functions on the sphere as the representation spaces. In the notation of 
Gelfand, Graev and Vilenkin, representations are labeled by the two numbers (m, n 2 ) or by (s, w), 
with (m, n 2 ) = {w — s + 1, w + s + 1). The 's' is the same 's' as in the spin weighted functions and 
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l w' is the conformal weight [33] (sometimes called 'boost weight'). The different representations 
are written as Dr niin2 y The special case of irreducible unitary representations, which occur when 
(ni,n 2 ) are not integers, play no role for us and will not be discussed. We consider only the case 
when (711,712) are integers so that the (s,w) take integer or half integer values. If ri\ and n% are 
both positive integers or both negative integers, we have, respectively, the positive or negative 
integer representations. The representation space, for each (s,w), are the functions on the sphere, 
V(s,w)(Ci C)> that can be expanded in spin-weighted spherical harmonics, s Yi m (£, £), so that 

oo 

V(s,w) (C, = ^2v(lm) ,^m(C)C)- (340) 
l=s 

Under the action of the Lorentz group, Eq. (|335[) . they transform as 

V(s, w )(lC) = e lsX K w V ^ w) (C,C)- (341) 

These representations, in general, are neither irreducible nor totally reducible. For us the important 
point is that many of these representations possess an invariant finite-dimensional subspace which 
(often) corresponds to the usual finite dimensional tensor representation space. Under the transfor- 
mation, Eq. (|34ip . the finite number of coefficients in these subspaces transform among themselves. 
It is this fact which we heavily utilize. More specifically we have two related situations: (1) when 
the {n\,n2) are both positive integers, (or w > |s|), there will be finite dimensional invariant sub- 
spaces, Df >, which are spanned by the basis vectors s Y; m (C, C)i with I given in the range, |s| < I 
< w. All the finite dimensional representations can be obtained in this manner. And (2) when 
the (—Tii,— Ha) are both negative integers (i.e., we have a negative integer representation) there 
will be an infinite dimensional invariant subspace, D7_ n ^ _ n \, described elsewhere |21j . One, how- 
ever, can obtain a finite dimensional representation for each negative integer case by the following 
construction: One forms the factor space, D(_ ni! _ n2 )/£)7 ni _ n2 y This space is isomorphic to one 
of the finite dimensional spaces associated with the positive integers. The explicit form of the 
isomorphism, which is not needed here, is given in Held et al[21| [16], 

Of major interest for us is not so-much the invariant subspaces but instead their interactions 
with their compliments (the full vector space modulo the invariant subspace). Under the action 
of the Lorentz transformations applied to a general vector in the representation space, the com- 
ponents of the invariant subspaces remain in the invariant subspace but in addition components 
of the complement move into the invariant subspace. On the other hand, the components of the 
invariant subspaces do not move into the complement subspace: the transformed components of 
the compliment involve only the original compliment components. The transformation thus has a 
non-trival Jordan form. 

Rather than give the full description of these invariant subspaces we confine ourselves to the 
few cases of relevance to us. 

I. Though our interest is primarily in the negative integer representations, we first address the 
positive integer case of the s — and w = 1, [(tii,^) = (2,2)], representation. The harmonics, 
I = (0, 1) form the invariant subspace. The cut function, X(r, £, £), for each fix values of r, lies in 
this space. 
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We write the GCF as 

u = X(C(t),C,C), (342) 

l,\m\<l 

After the Lorentz transformation, the geometric slicings have not changed but their description 
in terms of (u, () has changed to that of (u',(',( ). This leads to 

u' = KX, (343) 
= X', 

= ex(cX)+ E H /im (e a )Y lm (cX). 

l,\m\<l 

Using the transformation properties of the invariant subspace and its compliment we see that 
the coordinate transformation must have the form; 

£/a =F a (£ b ,H lm (t b ),...), (344) 

in other words it moves the higher harmonic coefficients down to the I — (0, 1) coefficients. The 
higher harmonic coefficients transform among themselves; 

H' lm (£' a ) = F lm (...,H l ' m '(£ b ),...). (345) 
Treating the £ a and £' a as functions of r, we have 

V = d T X = v a X, a = v a l a + v a J2 Hlm ^ (£°)*WC,C) , (346) 

i,|m|<( 

ar ar ar 



V = KV, 

= «'° K+ E H' lm , a Y lm (C'X) } , 

\ l,\m\<l 

= v b F a , b ll' a + E H' lm , a Y lm (C,C) I • 

\ l,\m\<l J 

Our Jo-space coordinates, z a = £ a , and their r-derivatives, v a , are the coefficients of the I = 
(0, 1) harmonics in respectively the X and the V expansions. A Lorentz transformation induces a 
specific coordinate transformation and associated vector transformation on these coefficients. 

II. A second important example concerns the mass aspect, (where we have introduced the for 
simplicity of treatment of numerical factors) 

* = *(o,-3) = *° r o° + + V l3 Y 2 % + ... (347) 
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is an s — and w = —3, [(711,7*2) = (—2, —2)] quantity. The invariant subspace has a basis set of 
the harmonics with I > 2. The factor space is isomorphic to the finite dimensional positive integer 
space, [(tii, n 2 ) — (2, 2)] and hence the harmonic coefficients of I = (0, f ) lie in this finite dimensional 
representation space. From this isomorphism we know the function of the form, ^°Yg + ^YQ, 
have four coefficients proportional to the Bondi four-momentum 

P a = {Mc 1 P i ) 1 (348) 

a Lorentzian four- vector. Note that we have rescaled the ^° in Eq. (|347p by the Yq, differing from 
that of Eq. (|63j) in order to give the spherical harmonic coefficients immediate physical meaning 
without the use of the factors in equations Eq. (f64|) and Eq. ([65| . 

This is the justification for calling the / = (0, 1) harmonics of the mass aspect a Lorentzian 
four-vector. Technically, the Bondi four-momentum is a co-vector but we have allowed ourselves a 
slight notational irregularity. 

III. The Weyl tensor component, if)®, has s = 1 and w = —3, [(7*1,712) = (—3,-1)]. The 
associated finite dimensional factor space is isomorphic to the finite part of the s = —l,w — 
1, [(7*1,7*2) = (3, 1)] representation. We have that 

V>? = 4i,_ 3 ) = ^uYu + ^Y 2 % + ... (349) 
leads to the finite-dimensional representation space 

finite space = V'li^iT 1 - (350) 

The question of what finite tensor transformation this corresponds to is slightly more complicated 
than that of the previous examples of Lorentzian vectors. In fact it corresponds to the Lorentz 
transformations applied to (complex) self-dual antisymmetric two-index tensors |29j . We clarify 
this with an example from Maxwell theory: from a given E and B, the Maxwell tensor, F ab , and 
then its self-dual version can be constructed: 

\yab+ _ pab _|_ ^p*ab _ 

A Lorentz transformation applied to the tensor, W ab+ , is equivalent [30] to the same transformation 
applied to 

ipv* o(E + *B).. (351) 

These observations allow us to assign Lorentzian invariant physical meaning to our identifications 
of the Bondi momentum, P a and the complex mass dipole moment and angular momentum vector, 

^MassDipole + ^ • 

IV. Our last example is a general discussion of how to construct Lorentzian invariants from the 
representation spaces. Though we will confine our remarks to just the cases of s = 0, it is easy 
to extend them to non-vanishing s by having the two functions have respectively spin-weight s and 
— s. 

Consider pairs of conformally weighted functions (s = 0), Wr w \, Gr_ x0 _2), with weights respec- 
tively, (w, —w — 2). They are considered to be in dual spaces. Our claim is that the integrals of 
the form 



I = / G { ^ 2) W {w) dn (352) 
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are Lorentz invariants. 

We first point out that under the fractional linear transformation, £' -H> (, Eq([337|), the area 
element on the sphere 

dn = 4i * A 5 (353) 

(i + CO 2 



transforms as 



dQ' = K 2 dn. (354) 



This leads immediately to 

I = I G[_ w _ 2) W[ w) dn' (355) 
K- w - 2 G { ^ 2) K w W (w) K 2 dn, 
G { _ w _2)W {w) dtt, 
the claimed result. 

There are several immediate simple applications of Eq. p52[ ). By choosing an arbitrary w — 
—2, s — function, say G(_2)(Ci C)> we immediately have a Lorentzian scalar, 



N = J G ( - 2) dn = J G[_ 2) dCl'. (356) 

If this is made more specific by chosing G(_2) = V~ 2 , we have the remarkable result (proved in 
Appendix [DJ that this scalar yields the "H-space metric via; 

8n{g ab v a v b )- 1 = [ V~ 2 d£l. (357) 



A simple variant of this arises by taking the derivative of (|357[) with respect to v a , and multi- 
plying by an arbitrary vector, w a , leading to 

8i:w a v a (g cd v c v d )~ 2 = J V~ 3 w a Z, a dfl. (358) 
Many other versions can easily be found. 
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8 Discussion / Conclusion 

8.1 History /background 

The work reported in this document has had a very long gestation period. It began in 1965 |37) 
with the publication of a paper where a complex coordinate transformation was performed on 
the Schwarzschild/Reissner-Nordstrom solutions. This, in a precise sense, moved the 'point source' 
onto a complex world line in a complexified spacetime. It thereby led to a derivation of the spinning 
and the charged-spinning particle metrics. How and why this procedure worked was considered to 
be rather mysterious and a great deal of effort by a variety of people went into trying to unravel 
it. In the end, the use of the complex coordinate transformation for the derivation of these metrics 
appeared as if it was simply an accident, i.e., a trick with no immediate significance. Nevertheless, 
the idea of a complex world line, appearing in a natural manner, was an intriguing thought, which 
frequently returned. Some years later, working on an apparently unrelated subject, we studied and 
found the condition [6] for an NGC, in asymptotically-flat spacetime, to have a vanishing asymptotic 
shear. This condition (our previously discussed shear- free condition, Equation (|174p ). was closely 
related to Penrose's asymptotic twistor theory. In the flat-space case it lead to the Kerr theorem 
and totally shear-free NGCs. From a different point of view, searching for shear-free complex null 
surfaces, the good-cut equation was found with its four-complex parameter solution space. This 
lead to the theory of 'H-space. 

Years later, the different strands came together. The shear-free condition was found to be 
closely related to the good-cut equation; namely, that one equation could be transformed into the 
other. The major surprise came when we discovered that the regular solutions of either equation 
were generated by complex world lines in an auxiliary Minkowski space |26j . (These complex world 
lines could be thought of as being complex analytic curves in the associated "H-space. The deeper 
meaning of this remains a question still to be resolved, i.e., what role, if any, do these complex lines 
play in "H-space. So far we find none.) 

The complex world line mentioned above, associated with the spinning, charged and uncharged 
particle metrics, now can be seen as just a special case of these regular solutions. Since these metrics 
were algebraically special, among the many possible asymptotically shear-free NGCs there was (at 
least) one totally shear-free (rather than asymptotically shear-free) congruence. This was the one 
we first discovered in 1965, that became the complex center-of-mass world line (which coincided 
with the complex center of charge in the charged case.). This observation was the clue for how to 
search for the generalization of the special world line associated with algebraically-special metrics 
and thus, in general, how to look for the special world line (and congruence) to be identified with 
the complex center of mass. 

For the algebraically-special metrics, the null tetrad system at 3 + with one leg being the tangent 
null vector to the shear-free congruence leads to the vanishing of the asymptotic Weyl tensor 
component, i.e., ipQ — ipl — 0. For the general case, no tetrad exists with that property but one 
can always find a null tetrad with one leg being tangent to the shear-free congruence so that the 
I = 1 harmonics of vanish. It is precisely that choice of tetrad that led to our definition of the 
complex center of mass. 
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8.2 Other choices for physical identification 

The question of whether our definition of the complex center of mass is the best possible definition, 
or even a reasonable one, is not easy to answer. We did try to establish a criteria for choosing such 
a definition: (i) it should predict already known physical laws or reasonable new laws, (ii) it should 
have a clear geometric foundation and a logical consistency and (iii) it should agree with special 
cases, mainly the algebraically-special metrics or analogies with flat-space Maxwell theory. We did 
try out several other possible choices [29] and found them all failing. This clearly does not rule 
out others that we did not think of, but at the present our choice appears to be both natural and 
effective in making contact with physical phenomena. 

8.3 Predictions 

Our equations of motion are simultaneously satisfactory and unsatisfactory: they yield the equations 
of motion for an isolated object with a great deal of internal structure (time-dependent multipoles 
with the emission of gravitational and electromagnetic radiation) in the form of Newton's second 
law. In addition, they contain a definition of angular momentum with an angular momentum flux. 
The dipole part of the angular momentum flux agrees with classical E&M theory. Unfortunately, 
there appears to be no way to study or describe interacting particles in this manner. 

However, there are some areas where these ideas might be tested, though the effects would 
be very small. For example, earlier we saw that there was a contribution to the Bondi mass (an 
addition to the Reissner-Nordstrom mass) from the quadrupole moment, 

Mb - M RN = AM = -\gc- 6 Re Q^iMZ ■ (359) 

There were predicted contributions to both the momentum and angular momentum flux from the 
gravitational and electromagnetic quadrupole radiation as well as new terms in the definition of the 
angular momentum, e.g., charge/magnetic-dipole coupling term 

9 9 

^ 2 -2 il * -2 n i/ 

3? c Vi = 3<? c d m- 

There are other unfamiliar terms that can be thought of as predictions of this theoretical construct. 
How to possibly measure them is not at all clear. 

8.4 Summary of results 

1. From the asymptotic Weyl and Maxwell tensors, with their transformation properties, we were 
able (via the asymptotically shear-free NGC) to obtain two complex world lines - a complex 
'center of mass' and 'complex center of charge' in the auxiliary "H-space. When 'viewed' 
from a Bondi coordinate and tetrad system, this led to an expression for the real center of 
mass of the gravitating system and a kinematic expression for the total angular momentum 
(including intrinsic spin and orbital angular momentum), as seen from null infinity. It is 
interesting to observe that the kinematical expressions for the classical linear momentum and 
angular momentum came directly from dynamical laws (Bianchi identities) on the evolution 
of the Weyl tensor. 
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2. From the real parts of one of the asymptotic Bianchi identities, Equation (|5fi)) . we found the 
standard kinematic expression for the Bondi linear momentum, P — M£' R + . . . with extra 

terms ffr^lj' and the Mathisson-Papapetrou spin coupling, among others. The imaginary 
part was the angular momentum conservation law with a very natural looking flux expression 
of the form: 

J k 1 = Flux Grav + Fluxg&M Q uad + Fluxg &M dipolc 

with 

J = spin + orbital + precession + varying magnetic dipole + quadrupole terms. 
The last flux term is identical to that calculated from classical electromagnetic theory. 

3. Using the kinematic expression for the Bondi momentum in a second Bianchi identity, Equa- 
tion (|57|) . we obtained a second-order ODE for the center of mass world line that could be iden- 
tified with Newton's second law with radiation reaction forces and recoil forces, MbCr = F. 

4. From Bondi's mass/energy loss theorem we obtained the correct energy flux from the electro- 
magnetic dipole and quadrupole radiation as well as the gravitational quadrupole radiation. 

5. From the specialized case where the two world lines coincide and the definitions of spin and 
magnetic moment, we obtained the Dirac gyromagnetic ratio, 3 = 2. In addition, we find the 
classical electrodynamic radiation reaction term with the correct numerical factors. In this 
case we have the identifications of the meaning of the complex position vector: £ l = £ l R 

£r = center-of-mass position 

S l ~ Mc£} = spin angular momentum 

D l E — qt; R = electric dipole moment 

D l M — q£\ — magnetic dipole moment 

8.5 A variety of issues and questions 

1. A particularly interesting issue raised by our equations is that of the run-away (unstable) 
behavior of the equations of motion for a charged particle (with or without an external field). 
We saw in Equation (|328l) that there was a driving term in the equation of motion depending 
on the electric dipole moment (or the real center of charge). This driving term was totally 
independent of the real center of mass and thus does not lead to the classical instability. 
However, if we restrict the complex center of charge to be the same as the complex center of 
mass (a severe, but very attractive restriction leading to g = 2), then the innocuous driving 
dipole term becomes the classical radiation reaction term - suggesting instability. (Note that 
in this coinciding case there was no model building - aside from the coinciding lines - and no 
mass renormalization.) 

A natural question then is: does this unstable behavior really remain? In other words, is it 
possible that the large number of extra terms in the gravitational radiation reaction or the 
momentum recoil force might stabilize the situation? Answering this question is extremely 
difficult. If the gravitational effects do not stabilize, then - at least in this special case - 
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the Einstein-Maxwell equations are unstable, since the run-away behavior would lead to an 
infinite amount of electromagnetic dipole energy loss. 

An alternative possible resolution to the classical run-away problem is simply to say that 
the classical electrodynamic model is wrong; and that one must treat the center of charge as 
different from the center of mass with its own dynamics. 

2. The interpretation and analysis of the complex analytic curves associated with the shear- 
free and asymptotically shear free null geodesic congruences naturally takes place in 'H-space. 
Though we get extraordinarily attractive physical results - almost all coinciding with standard 
physical understanding - It nevertheless is a total mystery as to what it means or what is the 
physical significance of this complex Ricci-fiat four-dimensional space. 

3. In our approximations, it was assumed that the complex world line yielded a cut of 3 + that 
were close to Bondi cuts. At the present we do not have any straightforward means of finding 
the world lines and their associated cuts of 3 + that are far from the Bondi cuts. 

4. As mentioned earlier, when the gravitational and electromagnetic world lines coincide we find 
the rather surprising result of the Dirac value for the gyromagnetic ratio. Unfortunately, 
though this appears to be a significant result, we do not have any deeper understanding of it. 
It was simply there for us to observe. 

5. Is it possible that the complex structures that we have been seeing and using are more than 
just a technical device to organize ideas, and that they have a deeper significance? One 
direction to explore this is via Penrose's twistor and asymptotic twistor theory. It is known 
that much of the material described here is closely related to twistor theory; an example is 
the fact that asymptotic shear-free NGCs are really a special case of the Kerr theorem, an 
important application of twistor theory (see Appendix [Aj . This connection is being further 
explored. 

6. With much of the kinematics and dynamics of ordinary classical mechanics sitting in our 
results, i.e., in classical GR, is it possible that ordinary particle quantization could play a role 
in understanding quantum gravity? Attempts along this line have been made [TH [7] but, so 
far, without much success. 

7. We reiterate that, a priori, there is no reason to suspect or believe that the world lines 
associated with shear-free congruences would allow the choice of a special congruence - and 
a special world line - to be singled out - and that furthermore it would be so connected with 
physical kinematics and dynamical laws. These results certainly greatly surprised and pleased 
us. 

8. As a final remark, we want to point out that there is an issue that we have ignored, namely, do 
the asymptotic solutions of the Einstein equations that we have discussed and used throughout 
this work really exist. By 'really exist' we mean the following: Are there, in sufficiently general 
circumstances, Cauchy surfaces with physically-given data, such that their evolution yields 
these asymptotic solutions? We have tacitly assumed throughout, with physical justification 
but no rigorous mathematical justification, that the full interior vacuum Einstein equations do 
lead to these asymptotic situations. However, there has been a great deal of deep and difficult 
analysis p~3l [101 [11] showing, in fact, that large classes of solutions to the Cauchy problem 
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with physically-relevant data do lead to the asymptotic behavior we have discussed. Recently 
there has been progress made on the same problem for the Einstein-Maxwell equations. 

9. An interesting issue, not yet explored but potentially important, is what can be said about 
the structure of "H-space where there are special points that are related to the real cuts of null 
infinity. 
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A Twistor Theory 

Throughout this review, the study of the asymptotic gravitational field has been at the heart 
of all our investigations. Here we make contact with Penrose's asymptotic twistor theory (see, 
e.g., [51[ l48j [22]). We give here a brief overview of asymptotic twistor theory and its connection 
to the good-cut equation and the study of asymptotically shear-free NGCs at 3 + . For a more in 
depth exposition of this connection, see [36] . 

Let M. be any asymptotically- flat spacetime manifold, with conformal future null infinity 3 + , 
coordinatized by (ub, C> C)- We can consider the complexification of 3 + , referred to as 3^, which is 
in turn coordinatized by (its, C): where now ub G C and C is different, but close to £. Assuming 
analytic asymptotic Bondi shear ct (ub,C, CX it can then be analytically continued to Jjt, i.e., we 
can consider ct°(mb, C, 0- 

We have seen in Section [4] that solutions to the good-cut equation 

9 2 G = fj°(G,C,C) (360) 

yield a four complex parameter family of solutions, given by 

u B = G(z a ;C,0- (361) 

In our prior discussions, we interpreted these solutions as defining a four (complex) parameter 
family of surfaces on 3^ corresponding to each choice of the parameters z a . 

In order to force agreement with the conventional description of Penrose's asymptotic twistor 
theory we must use the complex conjugate good-cut equation 

g 2 G = 7f°(G,C,C), (362) 
whose properties are identical to that of the good-cut equation. Its solutions, written as 

UB -G(z Q ;C,0, (363) 

define complex two-surfaces in 3^ for fixed ~z a . If, in addition to fixing the z a , we fix £ = Co € C, then 
Equation (|362|) becomes an ordinary second-order differential equation with solutions describing 
curves in (wb,C) space. Hence, each solution to this ODE is given by specifying initial conditions 
for G and d^G at some arbitrary initial point, C = Co- Note that it is not necessary that Co = Co 
on Dfjjt. However, we chose this initial point to be the complex conjugate of the constant Co ; i-e., we 
take G and its first C derivative at C = Co as the initial conditions. 

Then the initial conditions for Equation (|362[) can be written as [3 6) 

ubo = G(Co,Co), (364) 

_ fjr< 

L = gG(Co,Co) = Po^ r (Co,Co), 
<9Co 

with Pq = 1 + CoCo- Asymptotic projective twistor space, denoted PT, is the space of all curves in 
generated by initial condition triplets (ubo, Co> Lq) [51J : an asymptotic projective twistor is the 
curve corresponding to (ubo, Co, Lq). 
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A particular subspace of PT, called null asymptotic projective twistor space (P9T), is the family 
of curves generated by initial conditions, which lie on (real) 3 + ; that is, at the initial point, Co = 
Co, the curve should cross the real U + , i.e., should be real, ug = ubo- Equivalently, an element of 
P9t can be said to intersect its dual curve (the solution generated by the complex conjugate initial 
conditions) at Co = Co- The effect of this is to reduce the three-dimensional complex twistor space 
to five real dimensions. 

In standard notation, asymptotic projective twistors are defined in terms of their three complex 
twistor coordinates, (fjP,^ 1 ^) |51) . These twistor coordinates may be re-expressed in terms of the 
asymptotic twistor curves by 

M° = u bq - ioCo, (365) 
M 1 = L + ( u B0 , 

C = Co- 

By only considering the twistor initial conditions Co = Co, we can drop the initial value notation, 
and just let ubo = «b and C = C- 

The connection of twistor theory with shear-free NGCs takes the form of the flat-space Kerr 
theorem [5T1 [36] : 

Theorem. Any analytic function on PT (projective twistor space) generates a shear-free NGC in 
Minkowski space. 

Any analytic function on projective twistor space generates a shear-free NGC in Minkowski 
space, i.e., from F(fi° , , C) = F(«b — L(,L + (ub,() = 0, one can construct a shear-free NGC 
in Minkowski space. The L — L(ub, C, C), which defines the congruence, is obtained by solving the 
algebraic equation 

F(ub-LC,L + (ub,0 =0. 
It automatically satisfies the complex conjugate shear-free condition 

3I + H = 0. 

We are interested in a version of the Kerr theorem that yields the regular asymptotically 
shear- free NGCs. Starting with the general four-parameter solution to Equation (|362[) . i.e., ub = 
G(z a ; C, C), we chose an arbitrary world line z a = C a (r), so that we have 

u B = G(C>),C,C)=*(t,C,C), (366) 
L(t, C, C) = 3( T )X(r,C,C)- 
By inserting these into the twistor coordinates, Equation (|365l) . we find 

M °(r,C,C) = u B -L( = X-{B (T) X, (367) 
^(r,C,C) = L + Cub = d {T) X + (X. (368) 

The [iP and fi 1 are now functions of r and C: the C is now to be treated as a fixed quantity, the 
complex conjugate of C, and not as an independent variable. 

By eliminating r in Equations (j36T[) and p68l) . we obtain a single function of \i , and 
C: namely, F(p°, fj, , C) = 0. Thus, the regular asymptotically shear-free NGCs are described by 
a special class of twistor functions. This is a special case of a generalized version of the Kerr 
theorem [5H I31)] . 
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B CR Structures 

A CR structure on a real three manifold Af , with local coordinates x a , is given intrinsically by 
equivalence classes of one-forms, one real, one complex and its complex conjugate [31] . If we 
denote the real one-form by I and the complex one-form by m, then these are defined up to the 
transformations: 

I a(x a )l, (369) 

to -» f(x a )m + g(x a )l. 

The (a, /, g) are functions on Af: a is nonvanishing and real, / and g are complex function with / 
nonvanishing. We further require that there be a three-fold linear-independence relation between 
these one-forms [3Tj : 

I A to A to ^ 0. (370) 

Any three-manifold with a CR structure is referred to as a three-dimensional CR manifold. 
There are special classes (referred to as embeddable) of three-dimensional CR manifolds that can 
be directly embedded into C 2 . 

We show how the choice of any specific asymptotically shear-free NGC induces a CR structure 
on 3 + . Though there are several ways of arriving at this CR structure, the simplest way is to look 
at the asymptotic null tetrad system associated with the asymptotically shear-free NGC, i.e., look 
at the (l* a , m* a , m* a , n* a ) of Equation (|274j) . The associated dual one-forms, restricted to 3 + 
(after a conformal rescaling of to), become (with a slight notational dishonesty), 

I* = du B —=d( — -= d(, (371) 

l + CC l + CC 

d( d( 



l + CC' " l + CC' 

with L — L(ub, Ci 0? satisfying the shear-free condition. (This same result could have been obtained 
by manipulating the exterior derivatives of the twistor coordinates, Equation (|365[1 .) 
The dual vectors - also describing the CR structure - are 

M - P h + L lk^ +L lk' (372) 

d-d - d 

= P—+L- — =9 (u ,+L 



d 



Out 



Therefore, for the situation discussed here, where we have singled out a unique asymptotically 
shear-free NGC and associated complex world line, we have a uniquely chosen CR structure induced 
on 3+. 

To see how our three manifold, 3 + , can be imbedded into C 2 we introduce the CR equation [32] 

MK = (5 {ub) K + L^-K = 
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and seek two independent (complex) solutions, K\ = K 1 (ub,CX)^2 — K 2 (uq, £, £) that define 
the embedding of 3 + into C 2 with coordinates (K\,K2). 

We have immediately that K\ — C = x — iy is a solution. The second solution is also easily 
found; we see directly from Equation f)l T5|) [38 , 

<5 {Ub) T + LT = 0, (373) 

that 

T = T(u B ,(,(), 

the inverse to ub = X(tX,0i i s a CR function and that we can consider 3 + to be embedded in 
the C 2 of (t,C). 
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C Tensorial Spin-s Spherical Harmonics 

Some time ago, the generalization of ordinary spherical harmonics Yi m (£, Q to spin- weighted func- 
tions (s)Yim((, C) ( e -g-: [21} [TTl [40] ) was developed to allow for harmonic expansions of spin- weighted 
functions on the sphere. In this paper we have instead used the tensorial form of these spin- weighted 
harmonics, the tensorial spin-s spherical harmonics, which are formed by taking appropriate linear 
combinations of the ( s )Yi m ((,Q |43j : 

V s — \ ~* K sm V, 
1 li...k — / J jx li...k(s) 1 lm, 

where the indices obey \s\ < Z, and the number of spatial indices (i.e., i...k) is equal to I. Explic- 
itly, these tensorial spin-weighted harmonics can be constructed directly from the parameterized 
Lorentzian null tetrad, Equation (|74|) : 



i a = 2(1 "f cc - } (* + CC, C + C, »C-»C,-i + CC) , (374) 
v " (l + CC,-(C + C),*C-iC,i + CC). 



m 



2(1 + CC) 
V2 



(o,i-c 2 ,-i(i + C 2 ),2C) 



2(1 + CC) 

p = l + CC- 

Taking the spatial parts of their duals, we obtain the one-forms 



^(c+c,-*(c-c),-i + cc). ( 375 ) 
-^(c+c, -*(c+a-i+cc), 

_1 ^i-C 2 ,-i(i + C 2 ),2C), 



a/2P 

?i - n>i = -V2ieijkmjfhk 



From this we define k as 



Y i\...k = m l m J ...m k , (376) 

Y u...k = fhifhj...fhk- 

The other harmonics are determined by the action of the 3-operator on the forms, Equation (|375[) . 
(with complex conjugates) via 

& = m, (377) 

3rn = 0, 

oVi = — ra, 

9c = 2m, 

3to = n — I = — c. 
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Specifically, the spin-s harmonics are defined by 

= §'" s (YL.k) , (378) 
= ^ (^H...fc) ■ 

We now present a table of the tensorial spherical harmonics up to I = 2, in terms of the tetrad. 
Higher harmonics can be found in |43j . 



1 = 

Y° 
r o 




= 1 


I = 


1 




= mi, 




= ~Ci 




1 = rhi 



I = 2 

Y£j = m l m j , Y^j = - [c l m 2 + m l c j ), = 3c iCj - 2<% 

Y 2ij = m t rhj, Y^ij = - (cifhj + rh lCj ) 

In addition, it is useful to give the explicit relations between these different harmonics in terms 
of the 9-operator and its conjugate. Indeed, we can see generally that applying <3 once raises the 
spin index by one, and applying 8 lowers the index by one. This in turn means that 

Other relations for I < 2 are given by 

3*2 = -2Y,\, 

W° = -2Y-\ 



ciY 2 — Y 1 

B2y2 _ yO _ ?)2y-2 
u I 2ij — I 2ij — u I 2ij ' 

^Y 2i j = —GY^j, 

Finally, due to the nonlinearity of the theory, we have been forced throughout this review to 
consider products of the tensorial spin-s spherical harmonics while expanding nonlinear expressions. 
These products can be expanded as a linear combination of individual harmonics using Clebsch- 
Gordon expansions. The explicit expansions for products of harmonics with I = 1 or I = 2 are 
given below (we omit higher products due to the complexity of the expansion expressions). Further 
products can be found in [43l [29] . 



TENSORIAL SPINS SPHERICAL HARMONICS 
.1 Clebsch— Gordon expansions 



83 



I = 1 with I = 1 



Y 1 Y° - 


7f 


i 1 i 

ijkY lk + 2*2ij 


yly-1 _ 






yOyO _ 


5% 


+ -Y° 



I = 1 with I = 2 

yly2 t^3 

1 li I 2ij — J 3ijfc; 

^li ^2 j fc = — g ^ ( ^li ^2 j fe ) ' 



Z = 2 with Z = 2 

The Clebsch-Gordon expansions involving two I — 2 harmonics have been used in the text. 
They are fairly long and are not given here but can be found in [43]. 
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D %-space metric 

In the following the derivation of the H-space metric, is given. 

We begin with the cut function, u — Z (£ a (t) X , C) — X(t,(,C) that satishes the good cut 
eqation B 2 Z = a(Z, C, £)• The (£, £) are (for the time being) completely independent of each other 
though £ is to be treated as being "close" the complex conjugate of £. 

(Later we will introduce C* instead of £ via, 

c + w 

1-WQ 

for the purpose of simplifying an integration.) 

Taking the gradient of Z(z a , £, (), multiplied by an arbitrary four vector v a , (i.e., V = v a Z, a ), we 
see that it satisfies the linear Good Cut equation, 

Q 2 Z, a = a, z Z, a (380) 

s 2 y = <j,z v. 

Let Vo be a particular solution, and assume for the moment that the general solution can be written 
as 

Z, a = V Q l* a (381) 

with the four components of I* to be determined. Substituting Eq. (|38ip into the linearized GCE 
equation we have 



3 2 (^oC) 


= &,z 


VI* 


B(l* a B(V ) + V ®1) 


= <?,z 


vi* a , 


i:d 2 (Vo) + 2W m* a + v d 2 i* a 


= &,z 


V l* a , 


2dv m* a + v q% 


= o, 




2v w m* a + v 2 Q 2 r a 


- o, 




dv 2 m:+v 2 <5x 


= o, 




3(v 2 m:) 


= o, 





which integrates immediately to 

V 2 m* a = m* a (382) 

The m* are 3 independent 1 = 1, s = 1 functions. By taking linear combinations they can be 
written as 

m* a = T b a m b = T b a m b 

where t a is our usual l a = ^ ^1, — f^^, ~ ; iTcf ) ' ^ e coemc i ents are functions only of 
the coordinates, z a . 

Assuming that the monople term in V 2 is sufficiently large so that it has no zeros and then by 
rescaling V we can write V~ 2 as a monopole plus higher harmonics in the form 

V - 2 = 1 + dW, 
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where W is a spin-wt s = — 1 quantity. From Eq. (|382j) . we obtain 

K = v " 2 < = (i + ^R 

9C = m* a + BWm* a , 
Bl* a = m: + B(Wm* a ), 
&* a = Tffi b + T*B(Wm b ), 
m* a = T b d(l b + Wm b ), 

which integrates to 

r a =T b J a + Wm a ). 
The general solution to the linearized GCE is thus 

Z, a = V i: = VoT b {l a + Wm a ), 
V = v a Z, a =V v a T b (l a + Wm a ). 

We now demonstrate that 

GmvV)" 1 = (Sir)- 1 J v~ 2 dn. 

an = u^^. 
(i + CO 2 

In the integral, (|385|) . we replace the independent variables ((, £) by 

after some algebraical manipulation we obtain 
and (surprisingly) 



(f„ + Wm a ) = L a = — 1 



2 V i + CC i + C*C i + C*C 

so that 

Inserting Eas.([587|). (1555)) and (pMil) into ([3"8"5]) we obtain 



(Sa^V)- 1 - (Sir)- 1 J (Vov a T b L* b )- 2 V 2 dn* 
= (Sir)' 1 J (v a T b a L* b )- 2 dn\ 
= (Sir)- 1 f (v* b L* b )- 2 dn*. 
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Using the form Eq. (|389[) the last integral can be easily evaluated (most easily done using 9 and <p) 
leading to 

{9abV a v h )- x = (vabv* a v* b y 1 = (T c a T^ cd v a v b )-\ (392) 

gab = T^T b d 7J cd , 

our sort for relationship. 

We can go a step further. By taking the derivative of Eq. p91j) with respect to v a , we easily find 
the covariant form of v, namely 

. Va ^ = , 9abVb IY? = (Sit)- 1 { (v a T b L* b )- 3 T b Lldn* . 
{gabv a v b ) 2 [g a bv a v b ) 2 K ' J K a b ' a b 
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